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ABSTRACT 

We present the spectral energy distributions (SEDs) of a hard X-ray selected sample. The 
sample contains 136 sources with F2~wkeV >10~ 14 erg cm -2 s _1 and 132 are AGNs. The sources 
are detected in a 1 deg 2 area of the XMM-A'ewiora-Medium Deep Survey where optical data from 
the VVDS, CFHTLS surveys, and infrared data from the SWIRE survey are available. Based 
on a SED fitting technique we derive photometric redshifts with ct(1+z)=0.11 and 6% of outliers 
and identify AGN signatures in 83% of the objects. This fraction is higher than derived when a 
spectroscopic classification is available. The remaining 17t@% of AGNs shows star- forming galaxy 
SEDs (SF class). The sources with AGN signatures are divided in two classes, AGN1 (33^%) 
and AGN2 (50±?j%). The AGN1 and AGN2 classes include sources whose SEDs are fitted by 
type 1 and type 2 AGN templates, respectively. On average, AGNls show soft X-ray spectra, 
consistent with being unabsorbed, while AGN2s and SFs show hard X-ray spectra, consistent 
with being absorbed. The analysis of the average SEDs as a function of X-ray luminosity shows a 
reddening of the IR SEDs, consistent with a decreasing contribution from the host galaxy at higher 
luminosities. The AGNs in the SF classes are likely obscured in the mid-infrared, as suggested 
by their low £3-20^™/ '^o°5-iokeV ra tios. We confirm the previously found correlation for AGNs 
between the radio luminosity and the X-ray and the mid-infrared luminosities. The X-ray-radio 
correlation can be used to identify heavily absorbed AGNs. However, the estimated radio fluxes 
for the missing AGN population responsible for the bulk of the background at E >10 keV are 
too faint to be detected even in the deepest current radio surveys. 

Subject headings: X-rays: galaxies - Infrared: galaxies - quasars: general - Galaxies: active 
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1. Introduction 

AGNs are mostly found in X-ray surveys where 
they represent the majority of the detected pop- 
ulation. However, some AGNs are missed even in 
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the deepest X-ray observati ons, e.g. the Chan- 
dra Deep Surveys , CDFN dBrandt et all 12001 
Barger et al.l l2002|) and CPFS (jGiacconi et al 
2002: iRosati et all l200aY The existence of a 



large population of AGNs that is still undetected, 
even in the deepest X-ray surveys, i s suggested by 



the u nresolved X-ray background (jWorsley et al 



2005) and by the discovery of AGNs not de- 
tected at X-ray energies, b ut identified at radio 
and infrared wave l engths dPolletta et al 200G; 
Higdon et all 120051: iDonlev et al.l 120051) . Multi- 
wavelength observations of AGNs are necessary to 
make a census of AGNs and to fully characterize 
them. 

Current facilities and projects permit to study 
large AGN samples throughout the entire wave- 
length spectrum, from radio wavelengths to X-ray 
energies. Several wide and deep multi-wavelength 
surveys are currently being performed providing 
enough area, coverage and sensitivity to sample 
the mult i- wavelength properties of AGNs. A dif- 
ficult task that multi-wavelength surveys need to 
face is the identification of AGN. AGNs usually 
represent only a small fraction of sources com- 
pared to the population detected at optical and 
IR wavelengths and their properties can be elu- 
sive. Obscured AGNs can be elusive because their 
optical emission resembles that of normal galax- 
ies and therefore they cannot be identified by 
the color-selection techniques used for QSOs (e.g. 
Schmidt & Green 1983), or because their spec- 
tra do not show any AGN signature. Their IR 
emission can be dominated by PAH features, typ- 
ical of star-forming galaxies, instead of continuum 
emission from hot dust. Radio observations re- 
veal the radio-loud population of obscured AGNs 
(e.g. radio galaxies), but this comprises only a 
small fraction, ~10%, of the total AGN popula- 
tion. A hard X-ray selection fails to detect the 
luminous narrow-line AGNs in the numbers pre- 
dicted by some models of the CXRB. 

In order to obtain a complete picture of the 
AGN population, selection biases need to be taken 
into account. The fraction of AGNs missed in ob- 
servations in narrow energy windows, e.g. X-rays, 
infrared, or optical, is not well known. Moreover, 
the dispersion in the spectral energy distributions 
(SEDs) of all kind of AGNs needs to be well char- 
acterized to evaluate which and what fraction of 
AGNs are missed in a wavelength window and de- 



tected in another because of the unmatched depths 
of the observations at different wavelengths. The 
SED s are also important for evolutionary mod- 
els dTreister et alj2004:ISilva et alj2004 IXu et al 



2003; iBallantvne et al.l I2006T ). and to estimate 



the contribution of accretion to the cosmic back- 
grounds. 

In this work, we characterize the SEDs of 
a sample of AGNs selected in the hard X-rays 
and investigate how they vary as a function of 
X-ray absorption, luminosity, redshift, and host 
galaxy contribution. A selection in the hard X- 
rays minimizes the bias against absorbed AGNs 
compared to samples selected at soft or broad- 
band X-ray energies. Our analysis is based on 
the sources detected by XMM-Newton in a 1 
deg 2 region of th e XMM-Newton Mediu m Deep 
Survey (XMDS; IChiappetti et all I2005L , here- 
inafter Paper I) in th e XMM-Large Sca l e Sur - 



vey field (XMM-LSS; iPierre et al.ll2004 120071) . 



This region has also been observed as part of 
the Spitzer Wid e- Area Infrared Extragalactic Sur- 
vey (SWIRE) (|Lonsdale et all l2003h legacy pro- 
gram, of the VIMOS V L T Deep Survey ( WPS ; 
McCracken et al.l 120031: iLe Fevre etail |2004 . 
and of the Canada-France Hawaii Telescope 
Legacy "Peep" (PI) Survey. The available multi- 
wavelength data set, the source catalog, and the 
X-ray observations and data are summarized in 
Section [2] and described in detail in a companion 
paper (jTaier et al.l 120071 hereinafter Paper II) . 
Here, we present the technique used to classify 
the SEPs and to estimate photometric redshifts 
(Section [3]). We divide the sample in three main 
AGN classes based on their SEPs: AGN1, AGN2, 
and star-forming (SF) like AGNs. The average 
properties of the SEPs of the three AGN classes 
are described in Section |4j Luminosities and ab- 
sorption are estimated and their role in shaping 
the observed SEPs are discussed in Section [5] 
The relationship between the X-ray and the mid- 
IR luminosities is investigated in Section [51 and 
the variations of the mid-IR over X-ray luminosity 
ratio as a function of X-ray luminosity and AGN 
class are discussed in Section[6] The average SEPs 
of the three classes regrouped by X-ray luminos- 
ity are discussed in Section [JJ The analysis of 
the AGN IR colors compared with their luminos- 
ity, redshift and SEP classification is presented in 
Section [8] For a sub-set of sources, radio data are 
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also available. The known radio-IR and radio- X- 
ray correlations for AGNs and normal galaxies are 
compared with our data in Section [9] Our results 
are discussed in Section 1101 and summarized in 
Section [HQ 

Throughout the paper, we adopt a flat cosmol- 
ogy with H n = 71 kms" 1 Mp c" 1 , S1 M =0.27 and 
Q. A =0.73 ijSpergel et alJl2003h . For brevity, we re- 



fer to the sources in the sample by their XMDS 
sequence number (XMDS seq). The names of 
the sources, complying the IAU standard, along 
with the associated identifiers, are reported in Ta- 
ble [T] We use the terms red and blue to de- 
scribe the broad-band SED shape, where red and 
blue mean, respectively, increasing and decreasing 
fluxes, in vF v , at longer wavelengths. Fluxes re- 
fer to observed, not absorption or k-corrected val- 
ues, unless specified otherwise. The term absorbed 
refers to X-ray sources with effective column densi- 
ties Nh>10 22 cm~ 2 , and obscured to sources with 
red optical and IR spectra implying dust extinc- 
tion. The X-ray luminosity refers to the rest-frame 
broad band (0.5-10 keV) luminosity corrected for 
intrinsic absorption, unless specified otherwise. 
We use the terms type 1 and type 2 AGNs for 
AGNs with broad and narrow optical emission 
lines, respectively, in their optical spectra. We 
introduce a new nomenclature for the broad-band 
SED classification: AGN1, AGN2, and SF. The 
term AGN1 refers to an SED that is fitted by a 
type 1 AGN template, the term AGN2 refers to an 
SED that is fitted by a type 2 AGNs or composite 
template, and the term SF, which stands for star- 
forming, refers to an SED that is fitted by a star- 
forming galaxy template, like spirals and starburst 
galaxies. Although the templates are classified on 
the basis of their optical spectral properties, we do 
not expect a perfect match between the SED clas- 
sification and the properties of the optical spectra. 
The SED classification is more indicative of the 
main emitting components dominating at optical 
and IR wavelengths rather than of their optical 
spectra. Moreover, optic al spectra can be aff ected 
by host g alaxy dilution ( Moran et al.l 2002), ob - 
scuration (jBarger et al.l l200ll : iRigbv et al.l 120061 ) , 
instrumental and observational limitations, e.g. 
the lack of emission lines for sources in certain 
redshift ranges. 



2. Sample and multi-wavelength data set 

The sample studied in this work includes 136 
sources detected at > 3<7 in the 2-10 keV X-ray 
band, corresponding to ~10~ 14 ergcm _2 s _1 , in 
a 1 deg 2 region of the XMDS. The XMDS con- 
sists of 19 XMM-Newton pointings of 20 ksec 
nominal exposure covering a contiguous area of 
about 2.6 deg 2 . The XMDS survey is complete 
to a 2-10 keV flux of 9xl0" 14 ergcm- 2 s" 1 (Pa- 
per I). The 1 deg 2 area chosen for this study 
is centered at 02000 = 02 /l 26 m , and #2000 = 
— 04°30' (see Paper II for the layout of the se- 
lected field on the entire XMDS area) and bene- 
fits from a wealth of observations from various 
surveys: XM M-Newton data f r om th e XMM- 
LSS Survey (jPierre et all I2004L l2007t ): optical 



multi-band broad photometric data from the VI 



MOS VLT Deep Survey (V VDS; iMcGracken et al 



2003; ILe Fevre et al.ll2004f) , and from the Canada- 
France Hawaii Telescope Legacy "Wide" (Wl) and 
"Deep" (Dl) Surveys (CFHTLS^); near-infrared 
broad band photometric data from t he UKIRT 
Infrar ed Deep Sky Survey (UKIDSS; iDve et al 



2006; Law rence et al. 120061) : infrared data from 



the Spitzer Wi de- Area Infrared Extra galactic Sur- 
vey (SWIRE; lLonsdale et"aH l2003j); and radio 



data from the VLA VLM OS Survey (|Bondi et al 
2003; ICiliegi et al. 2005 ). Optical spectroscopic 



data are also available from the VVDS surve 



for a small sub-sample ( Le Fevre et al.l 1200 



5; 



Gavignaud et al.l |2006), from observations of a 



select ed sample of IR-selected AGNs (jLacv et al 



20061 ). and from the 2dF survey (Garcet et al., in 
prep.). 

The X-ray observations and data are described 
in Paper I, and the multi-wavelength catalog and 
data set are presented in Paper II. The X-ray data 
cover the energy range from 0.3 to 10 keV, and 
fluxes are measured in 5 different bands, 0.3-0.5, 
0.5-2, 2-4.5, 4.5-10, and 2-10 keV. Column den- 
sities are estimated from the hardness ratio HR 
derived from the soft (0.5-2 keV) and hard (2- 
10 keV) counts. For 55 sources with more than 
50 net counts, the column density is e stimated 
throu gh spectral modeling with XSPEC (jArnaud 
19961) (for more details see Paper II). 

The VVDS data set includes broad-band pho- 



1 http://www.cfht. hawaii.edu/Science/CFHTLS/ 
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tometric data in the B, V, R, and I niters, to 
an AB limiting magnitude (50% completeness for 
point s ources) of 26.5, 26 . 2, 25. 9, and 25.0, respec- 
tively (jMcCracken et al.l l2003h . Imaging data in 
the U band to an AB limiting magnitude of ~25.4 
are also available. The CFHTLS data set includes 
broad-band photometric data in the u* , g', r', i', 
and z filters, to an AB limiting magnitude (50% 
completeness) of i =24.5 in the wide survey and 
z'=26.1 in the deep survey. 

Near-infrared data in the J and K bands are 
available for a small fraction (~4%) of the area 
down to an AB limiting magnitude (50% com- 
pleteness) of 24.2 and 23.9. Additional near- 
infrared data from the UKIDSS project are also 
available. The UKIDSS observations cover about 
54% of the XMDS area and include J and K band 
data to a limiting magnitude of 22.3 and 20.8, re- 
spectively. 

The Spitzer infrared data include broad-band 
photometric data from the four IRAC bands at 
3.6, 4.5, 5.8, and 8.0 /im to a 5tr limit of 4.3, 8.3, 
58.5, and 65.7 /zJy, respectively, and MIPS data 
at 24 /jm to a 5a limit of 241 //Jy. The SWIRE 
observations cover about 90% of the area and are 
available for 122 sources, 91% of the entire sample. 

The radio data are from the VLA VIMOS 
Survey which reaches a 5a depth of 80 /zJy at 
1.4 GHz (jBondi et al.ll2003t ICiliegi et alJl2005h . 

The associations between the X-ray source list 
and the multi-wavelength catalogs are described 
in detail in Papers I and II. The selected multi- 
wavelength catalog contains 136 sources. Two 
sources are galaxy clusters (see Paper II) and 
are not included in the following analysis. Of 
the remaining 134 sources, 130 have optical coun- 
terparts, unique for 121 and multiple candidates 
for the other remaining 9 sources. Out of four 
sources without optical counterpart, one is close 
to a bright star and the photometry could not be 
extracted and the remaining 3 sources are opti- 
cally blank fields. All of the 122 sources that fall 
into the SWIRE area are detected in at least two 
Spitzer bands, and 87 are also detected at 24 /jm. 
UKIDSS data are available for 66 sources, of which 
55 are also SWIRE sources and 64 have an opti- 
cal counterpart. Radio data are available for 32 
sources, of which 30 are also SWIRE sources and 
30 have an optical counterpart. 



3. Template fitting &: photometric redshift 
technique 

In order to classify the spectral energy distribu- 
tions (SEDs) and estimate photometric redshifts 
of the sources in the sample, optical and IR data 
are combined and fitted with a library of galaxy 
and AGN templa tes. The SEDs are fit ted using 
the Hyper-z code ( Bolzonella et al. 20001 ). Hyper- 
z has been already successfully applied to vari- 
ous data sets and offers the possibility of using 
any template library and filters set. Hyper-z finds 
the best-fit by minimizing the \ 2 derived from 
the comparison of the observed SED and expected 
SEDs, at various redshifts, derived from the tem- 
plates using the same photometric system. The 
effects of dust extinction are taken into account 
by reddening the reference templates according to 
a selected reddening law. We use the prescrip- 
tion for e xtinction derived from high-redshift star- 
bursts bv lCalzetti et al. (2000) . In order to limit 
degeneracies in the best-fit solutions we limit the 
extinction Ay to be less than 0.55 mag and in- 
clude templates of highly extincted objects to fit 
heavily obscured sources. 

The library contains 20 templates including 1 
elliptical, 7 spirals, 3 starbursts, 6 AGNs, and 
3 composite (starburst+AGN) templates covering 
the wavelength range between lOOOA and 1000/im. 
The elliptical, spiral and starburst tem plates were 
gener ated with the GRASIL code (jSilva et al 



1998). The 7 spirals range from early to late types 
(SO-Sd), the starburst templates correspond to the 
SEDs of NGC 6090, M 82 and Arp 220. In all of 
the spirals and starburst templates the spectral 
region between 5 and 12 /im, where many broad 
emission and absorption features are observed, 
was replaced using observed IR spectra from the 
PHT-S spectrometer on the Infrared space Obser- 
vatory and from IRS on Spitzer. 

Templates of moderately luminous AGN, rep- 
resenting Seyfert 1.8 and Seyfert 2 galaxies, were 
obtained by combining models, broad-band pho- 
tometric data (NED), and ISO-PHT-S spectra 
(Schulz, private communication) of a random sam- 
ple of 28 Seyfert galaxies. The other four AGN 
templates include three templates representing 
optically-selected QSOs with different values of 
infrared/optical flux ratios (QSOl, TQSOl, and 
BQSOl) and one type 2 QSO (QS02). The 
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Fig. 1. — SEDs in vF v versus A of the template library used in this work. Elliptical, spiral and starburst 
galaxy templates are shown in the left panel and AGN templates in the right panel, the type 2 AGNs are 
the top 6 and the type 1 AGNs are the bottom 3. The names of the templates are annotated. 
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QS01 templates are derived by combining the 
SDSS quasar composite spectrum and rest-frame 
infrared data of a sample of 35 S PSS /SWIRE 



quasars ( Hatziminaoglou et al. 20051 ). After nor- 



malizing each SED in the optical, we derived three 
templates with the same optical spectrum but 
three different IR SEDs. The QSOl IR template 
was obtained from the average fluxes of all the 
measurements regrouped in wavelength bins, the 
TQSOl IR template was obtained from the high- 
est 25% measurements per bin, and the BQSOl 
IR template was obtained from the lowest 25% 
measurements per bin. The type 2 QSO tem- 
plate (QS02) was obtained by combining the ob- 
served optical/near-IR spe ctrum of the red qu asar 
FIRST J013435.7-093102 dGregg et alJl2002h and 
the rest-frame IR data from the quasars in the 
Palom ar-Green sample w ith consistent optical 
SEDs (|Polletta et al.ll2006l ). 



The composite (AGN+SB) templates are em- 
pirical templates created to fit the SEDs of the 
following objects: the heavily obscured BAL QSO 
Mrk 231 dBertall2005h . the Se yfert 2 galaxy IRA S 
19254-7245 South (119254; iBerta et alJ 120031 ). 
and the Sey fert 2 galaxy IRAS 22491-1808 (122491; 
Bert al 120051 ). All these objects contain a power- 
ful starburst component, mainly responsible for 
their large infrared luminosities (>10 12 L Q ), and 
an AGN component that contributes to the mid- 
IR luminosities. 

The full library of galaxy and AGN tem- 
plates is shown in Figure [TJ With respect to 
existing te mplate libraries deriv ed from empir- 
ical SEDs_|Colema^^L Il980h or from mod- 



els (iBruzual fc CharlqtJl2003l:lFioc fc Rocca-yolmerangel 
1997tlSilva et al.lll998HDevriendt et al.lll999l ). this 
library has a broader wavelength coverage and 
variety of spectral types. Examples of applica- 
tion of this library to va rious types of S WIRE 



sources can be fo u nd in lLonsdale et al.l (12 004'): 
iFranceschini et al I (120051): Hatziminaoglou et al] 
(l2005l):IJarrett et all (120061 ) : |Polletta et all (|2006l ): 
Weedman et al.l (|2006f) .~ 



Photometric rcdshift techniques have been 
applied to optical or optical and near-infrared 
data s ets providing reliable phot o metric red- 
shifts (fkowan- Robinson et al. 2005 : Ilbert et al 



20061 iBrodwin et alJ^OoiT However, with a lim- 



ited wavelength coverage, some spectral types can 
be degenerate, e.g., obscured starbursts and el- 



lipticals. In those cases, the inclusion of mid-IR 
data can break the degeneracy. However, the in- 
clusion of mid-IR data in photometric redshift 
techniques can degrade the photometric redshift 
estimates because the templates represent only 
a limited range of optical/IR ratios. We com- 
pared both approaches, fitting optical+near-IR 
and optical+near-IR-|- mid-IR data, and concluded 
that the latter method gives significantly better 
results both in the photometric redshift estimates 
and in spectral type classification. 

We fit the SEDs using the data from the op- 
tical to 24 ^m. Mid-IR (5.8-24 /mi) data are 
available for 107 out of the 134 sources in the 
sample. The best solution corresponds to the 
solution with minimum x 2 among all of the so- 
lutions with z=0-4 and B-band absolute mag- 
nitude within a pre-defined range. The allowed 
range of B-band absolute magnitude, Mb, varies 
with redshift and is different for normal galax- 
ies and for AGN templates. Mb must be greater 
than Mg m = — 6 — 5xLog(di), where dz, is 
the luminosity distance in Mpc, and lower than 
M maz = _o.7-5xLog(d L ), but always within the 
range —23.7, —17 for normal galaxy templates and 
greater than M™ = — 6.5 — 5xLog(d£) and lower 
than Mg" 1 = — 5xLog(di), but always within the 
range -28.8, -19 for AGN templates. The best-fit 
solution is the one with the lowest \ 2 among all of 
the solutions that satisfy the absolute-magnitude 
criterion. If there are no solutions in the defined 
absolute magnitude range, the criterion is modi- 
fied by extending the limits Mg"™ and M% ax by 
steps of 0.5 magnitudes, but always within the 
maximum allowed range, —23.7, —17 for normal 
galaxy solutions and —28.8, —19 for AGN solu- 
tions. The criterion was relaxed in 7 cases (see M 
flag in Table [J). 

In addition to the best solution, we consider up 
to two secondary solutions within the defined ab- 
solute magnitude limits. In most of the cases there 
are three solutions per source (the best and two 
secondary solutions), but in some cases, only one 
solution (or one minimum in the \ 2 distribution), 
is found within the defined absolute magnitude 
limits, thus for each source there might be from 
one to three solutions. We will refer to these as 
the acceptable solutions. The comparison between 
the spectroscopic and final photometric redshifts 
is shown in Figure [2l 
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Fig. 2. — Comparison between photometric 
and spectroscopic redshifts in Log(l+z) of the 
49 sources with available spectroscopic redshifts 
(open diamonds). Dashed lines represent 20% 
agreement in (1+z). The 5 catastrophic outliers 
are labeled with their XMDS seq: 48, 89, 439, 820, 
and 1259. Full circles (33 sources) represent the 
sources in the selected sample presented in this 
work. Crosses (16 sources) represent the addi- 
tional sources in the XMDS sample with spectro- 
scopic redshifts. Circles and plus signs are color- 
coded according to the SED class (blue: AGN1; 
red: AGN2; green: SF). A diamond is overplot- 
ted on each source and is color-coded according to 
the spectroscopic classification (blue: type 1 AGN; 
red: type 2 AGN; green: normal galaxy). The ver- 
tical dotted lines show the range of photometric 
redshifts of all acceptable solutions (note that the 
acceptable solutions can be at the extremes of the 
bar and not cover a continuous range in redshifts. 



The reliability and accuracy of the photometric 
redshifts are measured via the fractional error Az, 
the systematic mean error Az, the la dispersion 
a z , and the rate of catastrophic outliers, defined 
as the fraction of sources with \Az\ > 0.2. Az is 
defined as: 



Az = 



Zphot 



-spec 



and 



E 



Zphot 



l + z s 



/N 



(1) 



(2) 



with N being the number of sources with spectro- 
scopic redshifts. Since there are only 33 sources 
in the selected sample with a spectroscopic red- 
shift measurement (see Table [T]), we also included 
in the analysis of the photometric redshifts, all of 
the sources in the XMDS 4er VVDS sample (Paper 
I) with a spectroscopic redshift measurement, for 
a total of 49 sources. The systematic mean error, 
A~z, is -0.001 (-0.008), the rms, a z , is 0.12 (0.11), 
and the outlier fraction is 10% (6%) for the 49 
sources with spectroscopic redshifts. The values in 
parenthesis correspond to the results obtained us- 
ing only the 33 sources with spectroscopic redshifts 
in the hard X-ray selected sample. Although these 
results are not as satisfactory as those achieved 
for galaxy populations, even when few broad-band 



photo metric data are available (jBabbedge et al 
2004), they are better than what has been pre- 



viously obtained for AGN samples. The fraction 
of ou tliers for AGNs is indeed usually higher tha n 
25% (|Kitsionas et alJl2005t Eabbedge et al.ll2004 . 
The achieved accuracy does not allow us to per- 
form detailed analysis on single sources. However, 
it is adequate for a statistical analysis of the pop- 
ulation as presented in this work. Four (XMDS 
seq: 48, 89, 439, and 820) of the 5 outliers have 
| Az | < 0.28 and two of them (XMDS seq: 439, and 
820) are in the hard X-ray sample. The only out- 
lier with |Az| >0.28 is source XMDS seq 1259, 
which is not included in the selected sample. 

The redshift distributions of the whole sample, 
based on the photometric redshifts is shown with 
a dot-dashed line in Figure [3] The distribution 
peaks at z ~0.9, consistent with the reds hift distri- 
butio n of other X-ray selec t ed samples (Irlasinger 



20021: iBareer et al 
Eckart et alT l2006) . 
of the sources are 



20021 : ISzokolv et aiT l2004: 



About 54% (73 out of 134) 
at redshift below 1 and 90% 
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Fig. 3. — Redshift distribution of the entire sample 
(samples A+B), 134 sources (dash-dotted line), 
and of sample A only, 119 sources (dashed line) 
(see Section HJ. Each panel shows the redshift dis- 
tributions of the sources classified as AGN1 (top 
panel), AGN2 (middle panel), and SF (bottom 
panel) . The dotted lines correspond to the sources 
in each class in the entire sample, and the solid 
lines only to those in sample A. The shaded area 
represents the distribution of spectroscopic red- 
shifts in each class. 



are at redshift below 2.3. In the Figure, the red- 
shift distribution of the total sample is compared 
with that of AGN1, AGN2, and SF sources (see 
Section GTTJ). 

3.1. Spectral energy distributions 

According to the template corresponding to 
the best-fit solution selected by the procedure de- 
scribed above, each source is classified as one of 
the following broad classes, AGN1, AGN2, or SF. 
Examples of fits for the three classes are shown 
in Figure |H The AGN1 class includes all of the 
sources best-fitted with a QSOl template (bottom 
three on the right panel of Figure [TJ see also right 
panel of Figure [4]). They are characterized by blue 
optical SEDs (XF\ oc A~°' 9 '~ - 4 ) and red IR SEDs 
(XF\ oc A 2 ' 4 ). The change in slope, from blue to 
red, occurs at about 1 fmx. Emission from the ac- 
cretion disk dominates at shorter wavelengths and 
from AGN-heated dust at longer wavelength. The 
AGN2 class includes sources best-fitted with the 
Seyfert templates, the three composite AGN+SB 
templates, or the QS02 template (top six on the 
right panel of Figure [TJ and see also middle panel 
of Figure 0J). The AGN2 templates include a 
wide variety of SED shapes. They are mostly red 
at optical and near-IR wavelengths, but in some 
cases an upturn (blue) at short wavelengths (far- 
ultraviolet) is present. Their IR SEDs range from 
red power-law like SEDs, similar to the AGN1 or 
redder, to composite SEDs where emission fea- 
tures associated with star-formation are present 
in addition to continuum emission in the near-IR. 
The SF class includes sources fitted with a spiral 
or a starburst template (left panel of Figures [TJ 
and 2|) . SF templates are all red in the optical 
and blue in the near-IR, consistent with emis- 
sion from stellar populations. At mid-IR wave- 
lengths, the SED is dominated by emission fea- 
tures and a warm dust continuum associated with 
star- forming regions. Photometric and spectro- 
scopic redshifts, best-fit template and classifica- 
tion are reported in Table [JJ 

Although a best-fit template, and thus an SED 
classification, is available for all sources in the 
sample, in the following analysis we will discuss 
the results only for sources with a reliable fit and 
simply report the derived values for the rest of 
the sample for completeness. We define a fit re- 
liable if at least six data points could be fitted. 
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Fig. 4. — Scaled SEDs (full black circles and triangles) and best-fit templates (gray solid curves) in F v versus 
A of 18 sources: 6 classified as SF (left panel), AGN2 (middle panel), and AGN1 (right panel). Downward 
pointing arrows correspond to 5cr upper limits. The source sequence number and best-fit template names 
are listed in the same order as the SEDs are plotted. 
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Six sources (XMDS seq: 139, 199, 233, 403, 840, 
1219) have been observed only in the optical or in 
the IR and have less than six data points in total. 
Because of the limited wavelength coverage, they 
are not included in the following analysis. Nine 
sources with ambiguous counterpart (XMDS seq: 
91, 111, 120, 227, 255, 449, 470, 747, 844) are also 
excluded in the following analysis. We will refer to 
the 15 sources as sample B. In Table[T]we report all 
best-fit solutions for the sources with ambiguous 
counterpart in order of increasing distance of the 
optical counterpart, but in the following figures we 
will show only the best-fit solution corresponding 
to the closest counterpart. The sample that will 
be analyzed in the following sections contains the 
remaining 119 sources (the total 134 minus the 15 
in sample B; hereinafter sample A), of which 38 
are classified as AGNls (32% of the sample), 59 
(50%) as AGN2s and the remaining 22 (18%) as 
SFs. 

The reliability of this classification depends on 
the number of available optical and IR data points 
and on the ability of identifying the signatures that 
make the SEDs of the various classes unique. In 
the case of AGNls, the blue optical SED and red 
IR SED provide two clear signatures that make 
their classification highly reliable (see right panel 
in Figured]). The distinction between AGN2s and 
SFs is less obvious than the identification of AGN1 
since the only difference between their SEDs is in 
the shape of the near- and mid-IR SED (sec left 
and middle panels in Figure HJ. When only few IR 
data points are available or the data do not sample 
the region where the AGN signature is visible (i.g. 
3-5 /itm rest- frame), the separation between the 
two classes is more uncertain. In order to estimate 
the uncertainty of our classification we counted all 
sources with an acceptable solution in a different 
class (see Section [3] for a description of the accept- 
able solutions). For each class, we counted the 
number of sources with solutions in other classes, 
and the number of sources in other classes with 
a secondary solution in the specific class. For the 
AGN1 class, we estimate that 2 sources might be- 
long to a different class, AGN2 or SF, 4 AGN2s 
have an AGN1 secondary solution, and 2 SFs have 
an AGN1 secondary solution. Thus, there are 36 
sources with a secure AGN1 fit, and 44 sources 
with a possible AGN1 fit, which implies an AGN1 
fraction of 32±2%. Analogously, for the AGN2 



class, we estimate that 13 sources might belong 
to a different class, AGN1 or SF, 1 AGN1 has an 
AGN2 solution, and 5 SFs have an AGN2 solution. 
Thus, there are 46 sources with a secure AGN2 fit, 
and 65 sources with a possible AGN2 fit, which 
implies an AGN2 fraction of 50± 5 n %. For the SF 
class, we estimate that 7 sources might belong to a 
different class, AGN1 or AGN2, 1 AGN1 has an SF 
solution, and 9 AGN2s have an SF solution. Thus, 
there are 13 sources with a secure SF fit, and 30 
sources with a possible SF fit, which implies an 
SF fraction of 18+£%. Note that, by definition, 
the given uncertainties on the class fractions cor- 
respond to lower and upper limits to the fraction 
of sources in each class. 

The reliability and completeness of our classi- 
fication can be estimated using the spectroscopic 
classification, however we do not expect an exact 
correspondence between our SED-based classes 
and the spectroscopic classes. A spectroscopic 
classification is available for 49 sources, including 
sources that are not in the hard X-ray samples. Of 
these, 26 are spectroscopically classified as type 1 
AGNs, 12 are spectroscopically classified as type 
2 AGNs, and 11 as star-forming galaxies. Among 
the 26 type 1 AGNs, 16 are classified AGNls, 
9 AGN2s, and 1 SF, and all the AGNls with a 
spectroscopic classification (16) are type 1 AGNs. 
Thus, the SED classification is 62% complete (16 
AGNls out of 26 type 1 AGNs), and 100% reli- 
able (all AGNls are type 1 AGNs). Among the 
12 sources spectroscopically classified as type 2 
AGNs, 7 are classified as AGN2s, and the remain- 
ings as SFs. The AGN2 class (24 sources) includes 
9 type 1 AGNs, 7 type 2 AGNs, and 8 SFs. Thus, 
the completeness of the AGN2 sample is 58% (7 
AGN2 out of 12 type 2 AGNs), and the reliability 
is 29% (7 out of 24 AGN2s are type 2 AGNs). 
There are 11 sources with a spectroscopic clas- 
sification consistent with being normal galaxies. 
Eight of them are classified as AGN2s, and 5 as SF. 
Among the 9 sources classified as SF, 1 is a type 1 
AGN, 5 are type 2 AGNs, and 3 are star-forming 
galaxies according to the spectroscopic classifica- 
tion. The completeness of the SF sample is thus 
27% (3 out of 11 star- forming galaxies are classi- 
fied SF), and the reliability is 33% (3 out of 9 SF 
sources are star-forming galaxies). In summary, 
the SED and the spectroscopic classification agree 
in 26 of the 49 sources with spectroscopic classifi- 
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cation (53%). 

We can also compare the fraction of AGNs that 
are identified as AGNs by the SED classification 
(AGN1 or AGN2) and by the spectroscopic classi- 
fication (type 1 and type 2 AGNs). Note that all 
sources are AGNs as probed by their X-ray prop- 
erties. Out of the 49 AGNs with spectroscopic 
classification, 38 (78%) are spectroscopically iden- 
tified as AGNs (26 type 1 and 12 type 2), and 40 
(82%) are classified AGNs by the SED fitting tech- 
nique (25 AGN1 and 7 AGN2). The latter fraction 
is consistent with the 83% derived from the entire 
X-ray sample. 

3.2. Redshift-dependent derived quanti- 
ties: luminosities and absorption 

Using the available redshifts, spectroscopic if 
available, photometric otherwise, and the best-fit 
templates, we derive the luminosities in the broad 
and hard X-ray energy bands (0.5-10 keV and 2- 
10 keV rest-frame), in the mid-IR (3-20 /j,m rest- 
frame), and in the radio at 1.4 GHz rest-frame. 
The X-ray luminosities, hereinafter Lo.5_io keV 
and L2-10 keV i are derived assuming a power- 
law spectrum (F(E)oci? _ ^ r_1 ^) with photon in- 
dex r=2.0. Effective column densities are derived 
from the observed column densities after applying 
the redshift correction (see Paper II for details). 
In case the observed column density is consistent 
with the Galactic value, we derive an upper limit 
to the effective column by assuming the Galactic 
value and the redshift of the source (see also Sec- 
tion[5|). The observed column densities are used to 
correct the X-ray luminosity for absorption, here- 
inafter i§°5l 10 keV The uncertainty on the X-ray 
luminosities are estimated by propagation of the 
uncertainty on the flux and on the column density. 

The mid-IR luminosities, L^^oum, are esti- 
mated by integrating the best-fit template in the 
rest-frame 3-20 /xm wavelength range. To esti- 
mate the uncertainty on the mid-IR luminosities 
we first fit the IR data (X b s > 3.6/im) of each 
source with all of the templates, derive the x 2 and 
then use the range of luminosities obtained from 
all acceptable solutions. 

The radio luminosities, L1.4 GHz, are derived 
assuming a power-law spe ctrum (F„ oc v~ aR ) 
with spectral index au=0.8 (jCondon fc Broderick 
19881 ). The uncertainty on the radio luminosity is 



derived from the radio flux uncertainty. 

The X-ray, mid-IR and radio luminosities are 
listed in Table CD 

3.3. Star-forming galaxies 

In the following analysis, we will include only 
sources that contain an AGN. In order to identify 
sources whose X-ray emission is consistent with 
emission from a star- forming galaxy, e.g., from 
X-ray binaries, we require a broad X-ray lumi- 
nosity below 10 42 ergs _1 , and a soft X-ray spec- 
trum. Star-forming galaxies are typically charac- 
terized by luminosities below 10 erg s~ 1 and sho w 
steep X-ray spectral slopes dKim et al.l Il992al fbl: 
Colbert et all2004l:|Persic et al.ll2004h ~~This crite- 
rion might misidentify some normal star-forming 
galaxies with a large population of high-mass X- 
ray binaries ( HMXBs) that prod uces a flat X- 
ray spec trum dColbert et al. 20041) . (see also Fig- 
ure 2 in lAlexander et al.1 2005al ). and might also 
misidentify low-luminosity AGNs whose X-ray lu- 
minosity can be lower than 10 42 ergs _1 . 

There are 3 sources with an absorption-corrected 
luminosity below 10 42 ergs _1 , sources XMDS seq 
178, 842 (aka Arp 54), and 1248. Sources XMDS 
seq 178 and 1248 are characterized by low col- 
umn densities, 22xl0 20 cm -2 and 2.6xl0 20 cm~ 2 
(consistent with the Galactic value), and hardness 
ratios, —0.43, and —0.54, respectively. However, 
source XMDS seq 842 (z SJ , ec =0.043) has a hard- 
ness ratio of 0.16 which implies a column density 
of 1.78xl0 22 cm~ 2 . Because of its hard X-ray 
spectrum, source XMDS seq 842 is included in 
the AGN sample that will be investigated in the 
following sections. The classification of source 178 
is also dubious because its optical to X-ray flux 
ratio (Log(F . 5 _ 10 kev/Fji) ~ -0.49) is consis- 
tent with those typical of u nobscured AGNs (— 1< 
Log(F .^ w kev/F R ) < 1; lAkivama et al.ll2003l ). 
Indeed in Paper II, where star-forming galax- 
ies are identified based on the optical to X-ray 
flux ratio, XMDS seq 178 is considered an AGN 
and XMDS seq 842 is considered a star-forming 
galaxy. Here, we prefer to adopt a criterion which 
is based only on intrinsic X-ray properties (lumi- 
nosity and Njj) rather than flux ratios that are 
redshift-dependent. Thus, we consider XMDS seq 
178 and 1248 normal star-forming galaxies, and 
XMDS ID 842 an AGN. Sources XMDS seq 178 
and 1248 will not be taken into account in the 
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following analysis which will focus only on the 
AGN population (132 sources in total, and 117 in 
sample A). 

The fractions of AGNs in the three SED 
classes thus become 33+5_i% AGNls, 50+6_n% 
AGN2s, and 17+9_ 6 % SFs. 

4. Average SEDs of the three AGN SED 
classes 

After deriving the luminosity SEDs in vhv for 
each source assuming the source photometric red- 
shift, or spectroscopic redshift if available, aver- 
age normalized rest-frame SEDs, from the X-ray 
to the mid-IR, were obtained from all AGNs in 
the A sample (117 sources) divided by class (SF, 
AGN2 and AGN1). The average SEDs are shown 
in Figu re El Similarly to the median QSO tem- 
plate in I Elvis et al.l ( 19941 ). we normalize the SEDs 
at 1 /im rest-frame assuming the best-fit template 
for each source. This choice of wavelength is well 
suited to investigate separately the SED dispersion 
of the main emitting components in AGNs, the ac- 
cretion disk and the stars in the optical and the 
dust in the infrared. The average SEDs are derived 
from the normalized rest-frame SEDs by taking 
the mean of the logarithm of the wavelength and 
the weighed mean of the luminosity in Log(i/Li/) 
in 0.3 wide bins. Smaller bins were chosen if more 
than 40 data points were available or wider bins to 
include at least 10 points. The weights are based 
on the luminosity uncertainties. Upper limits in 
the optical and infrared are taken into account by 
assuming the upper limit to the luminosity and a 
logarithmic uncertainty equal to twice the maxi- 
mum logarithmic uncertainty of each source. Up- 
per limits to the X-ray luminosities are taken into 
account by assuming as luminosity the la lumi- 
nosity value and a logarithmic uncertainty equal 
to twice the maximum logarithmic uncertainty ob- 
served in the X-rays for that source. The average 
SF, AGN2, and AGN1 SEDs contain respectively 
15, 29, and 20 bins, and are obtained by averag- 
ing the rest-frame SEDs of 20 SFs, 59 AGN2s, and 
38 AGNls. The average SEDs are shown in Fig- 
ure [5] as thick black curves, 1 and 2a dispersion 
are shown as grey shaded areas. Since the redshift 
distribution of the sample is broad, up to z >3.5, 
and the majority of the sample has z the ob- 
served X-ray data sample rest-frame energies up 



Rest-frame Loq(i') (Hz) 
18 17 16 15 14 13 



^ 



-1 



-2 



2 - 



■i i i — ■ ■ ■ i ■ 



■i i ■ 



AGN1 (38) 

> > > ' ' 





AGN2 (59) 

I ' I ' ' 





SF (20) 



_i i i i_ 



0.0001 0.0010 0.0100 0.1000 1.0000 10.0000 
Rest-frame A (y^m) 



Fig. 5. — Normalized at 1 /im average rest-frame 
X-ray-mid-IR SEDs in \jOg{vLu) vs Log(A) (full 
circles connected by straight lines) of AGNs classi- 
fied as AGN1 (38 sources) (top panel), AGN2 (59 
sources) (middle panel) and SF (20 sources) (bot- 
tom panel). The light and dark shaded areas cor- 
respond to 2 and la dispersion, respectively. Note 
that the X-ray luminosities are not corrected for 
absorption. The black thin curves represent the 
median (dashed curve) ±90% disper s ion (d otted 
curves) QSO template bv lElvis et all (|l994l ). 
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to 45 keV for the highest- z sources, and up to 20 
keV for the majority of the sample. Thus, the av- 
erage SEDs sample rest-frame energies up to 20 
keV. 

The global SEDs, from X-ray to infrared, show 
a wide dispersion in each class. Since they are 
normalized at 1 /im, this dispersion is more pro- 
nounced at X-ray wavelengths. The average SEDs 
of the three classes show some clear differences. 
The average SEDs become increasingly blue in 
the optical-NIR (A rest < 1/zm), red in the IR 
(A res * ~ 1 — 10/im), and soft in the X-rays in the 
sequence SF^AGN2^AGN1 (from bottom to top 
in Figure[5|) . The average SEDs are compared with 
the med ian QSO templa te of optically-selected 
quasars (|Elvis et alJ Il994h in Figure \E\ The av- 
erage AGN1 SED is redder in the optical than 
Elvis's template, with only the bluest sources be- 
ing consistent with Elvis's template. The redder 
optical SED of the AGNls compared to Elvis's 
template suggests that these sourc es might in- 
clude some reddened type 1 AGNs (I Wilkes et al 
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Gregg et al. 
2005: 
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Glikman et al. 
Wilkes et al.ll2005l ). Indeed. 



Urrutia et al.l 
the observed range of extinction for AGNls in 
V-band is 0.0-0.55, and the median value is 0.40 
(see Table H}. The SF and AGN2 optical SEDs are 
almost identical and significantly redder than the 
AGN1 optical SEDs. The observed different shape 
suggests that the origin of the optical emission is 
different in AGN1 and in AGN2, and SF sources, 
rather than the effect of reddening on an AGN 
spectrum. The AGN2 and SF optical emission is 
likely dominated by the host galaxy emission in 
most of the cases. 

The near- and mid-IR (A ~l-20 /im) aver- 
age SEDs of AGNls and AGN2s are consistent 
with each other at la level, and both are redder 
than the SF average SED. The mid-IR SED of 
AGN1 and AGN2 are consistent with being domi- 
nated by AGN-heated hot dust emission, while, in 
SFs, the host galaxy dominates. The mid-IR SED 
of Elvis's template is consistent with the average 
AGN1 mid-IR SED, and redder than the average 
AGN2 mid-IR SED. 

In the X-rays, we show X-ray luminosities with 
no correction for intrinsic or Galactic absorption. 
AGNls show on average softer X-ray spectra than 
SFs and AGN2s, although there is a significant 
dispersion. The observed average spectra are con- 



sistent with increasing absorption in AGN2 com- 
pared to AGN1, and in SF compared to AGN2. 

Compared to Elvis's template, the normalized 
X-ray luminosity (Lx/L(l/xm)) of AGNls is more 
than two times lower. This difference might be 
due to a well known bias in Elvis's QSO sample 
which favors sources with strong X-ray-to-optical 
lumi nosity ratios and detections in the so ft X- 
rays (jWilkes k Elvidll987t lElvis et al.lll994l) . 

The observed differences in the optical and IR 
SEDs of the three classes are mainly a consequence 
of the spectral classification. However, the SED 
classification is independent of the X-ray proper- 
ties, therefore the differences in the X-ray average 
spectra indicate that harder X-ray spectra are ob- 
served in more obscured AGNs. 

5. Luminosities and X-ray absorption 
versus AGN SED class 

The absorption-corrected rest-frame broad- 
band X-ray luminosities are compared to the in- 
trinsic column densities in Figure [6] for the entire 
AGN sample (132 sources). Full symbols represent 
sources in sample A (117 sources), and open sym- 
bols represent sources in sample B (15 sources). 
Upper limits to the column density are adopted 
for sources with an observed column density lower 
or consistent with the Galactic value. The upper 
limit to the column density is derived from the 
Galactic value after applying the correction for 
redshift (N H Ga( x (1 + z) 2 - 6 ). Median luminosity 
values, derived from the A sample, for each class 
are shown with dashed lines. The absorption- 
corrected rest-frame X-ray luminosities of AGN1 
are on average more than ten times those of SF, 
while AGN2 show intermediate values. No corre- 
lation with the intrinsic absorption measured in 
the X-rays is observed. The paucity of heavily ab- 
sorbed (N H > 10 23 
below 10 44 ergs" 
because they would be too faint in the X-rays to 
be detected above our selection threshold. 

Note that the area (ldeg 2 ) of this survey pro- 
vides a large sample of highly luminous (Lg°5^ 10 keV 
10 45 ergs -1 ) AGNs (21 sources). As commonly 
found in X-ray selected samples where unabsorbed 
type 1 AGNs dominate the most luminous pop- 
ulation (jUeda et al.ll2003t lAkvlas et alj|2003 ). we 
also find a larger fraction of AGNls at higher X- 
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Fig. 6. — Rest-frame X-ray luminosities corrected 
for absorption as a function of the effective column 
density measured in the X-ray for the whole sam- 
ple (open and full circles). Full symbols represent 
sources in sample A and open symbols to sources 
in sample B (see Section [4]). Different symbols and 
colors correspond to the 3 SED classes, AGN1: 
blue circles, AGN2s: red reversed triangles, and 
SFs: green triangles. The dashed lines represent 
median values of the absorption-corrected X-ray 
luminosities for the 3 classes. 
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Fig. 7. — Comparison of mid-IR luminosity with 
the effective column density measured in the X- 
rays. Symbols as in Figure [6] 



ray luminosities. The fraction of AGNls increases 
from 22% at i^S-io keV < 1045 ergs^to 81% 
at L o°5-w keV > 1045 ergs" 1 . A detailed dis- 
cussion on the fraction of unobscured (AGNls) 
and obscured (AGN2s+SFs) AGNs as a function 
of X-ray luminosity and a comparison with pre- 
vious results is reported in Paper II. Paper II 
also discusses the uncertainty of these estimates 
due to the large corrections in luminosity and 
column density for sources at high-redshifts and 
the properties of absorbed QSOs (NH>10 22 cm -2 
and i^-To keV > 10 44 ergs -1 ). It is interesting 
to note the lack of SF sources in sample A at 
high X-ray luminosities, >10 45 ergs -1 . This in- 
dicates that powerful AGNs always dominate in 
the IR and are, therefore, identifiable as AGNs 
at those wavelengths. Thus, IR surveys can po- 
tentially be used to identify and study the most 
luminous AGNs. There are however two sources 
in sample B (XMDS seq: 403, 449) classified as SF 
at high luminosities (i^Tfo keV > l^ 45 er g s *)• 
Since they belong to the B sample, they are not in- 
cluded in the analysis (open triangles in Figure [6]). 
Deeper multi-wavelength observations are neces- 
sary to better characterize these sources (classi- 
fication and redshift) and verify their high X-ray 
luminosities. 

The mid-IR luminosity is compared with the 
column density measured in the X-rays in Fig- 
ure [3 Similarly to the distribution of the X-ray 
luminosities, the majority of AGNls show predom- 
inantly high mid-IR luminosities (>10 45 erg s" 1 ), 
while SFs show mostly low mid-IR luminosities 
(< 10 44 5 erg s" 1 ), and AGN2s show intermediate 
values (10 44 ~ 45 - 5 erg s" 1 ) overlapping with those of 
the other two classes. As for the X-ray luminosi- 
ties, the mid-IR luminosities of the objects in each 
class do not show any obvious correlation with the 
intrinsic absorption measured in the X-rays. 

In the next sections, we will investigate whether 
the observed multi-wavelength properties of the 
three classes can be simply explained by differ- 
ent luminosities, dust extinction or orientation ef- 
fects. These parameters will be investigated sepa- 
rately for each AGN class by analyzing the mid-IR 
and X-ray luminosity ratio, the dependency of the 
SEDs on the X-ray luminosity, and of the IR colors 
on the luminosity and dust extinction. 
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6. The range of mid-IR over X-ray lumi- 
nosity ratios 

The nuclear mid-IR continuum in AGNs is due 
to reprocessing of the AGN emission by circum- 
nuclear dust, e.g. the putative torus, and is thus 
function of both the AGN luminosity and the 
distribution of the obscuring matter. The dust 
covering factor in AGNs can be investigated by 
comparing the ratio between the mid-IR contin- 
uum (re-radiation) and the intrinsic X-ray emis- 
sion. A lower mid-IR/X-ray luminosity ratio is 
expected in highly luminous AGNs if the dust in- 
ner radius recedes (|Lawrencdll991l ) and the torus 
height is constant, and in type 2 sources where 
the torus is seen edge-on bec ause of mid-IR ab- 
sorption ( Granato et al.l fl997l ) . However, studies 
on Seyfert 1, Sey fert 2 galaxies ( Lutz et all 2004: 
Horst et al.ll2006[ ). and type 2 QSOs (|Sturm et al ' 



2006) where the dust covering factor is estimated 
from the L mid ^ IR /L^[ keV luminosity ratio, 
where L rn i ( i-iR is the monochromatic luminos- 
ity at 6 or 12 /im, have not found any dependency 
on AGN type, absorption or X-ray luminosity. 

In Figure [5J we compare the mid-IR/X-ray lu- 
minosity ratio with the absorption-corrected X- 
ray luminosity for the three AGN classes. In or- 
der to investigate whether there is a trend with 
the X-ray luminosity or the AGN class, we re- 
grouped the data for each class in five bins with 
the same number of sources per bin for each class, 
and derived the average ratio. The average val- 
ues are derived only for sample A (full symbols) 
and are shown as diamonds connected by solid 
lines. For clarity, we show the data of the single 
sources witho ut uncertainties. Analogously to pre- 
vious studies dLutz et all 12004 iHorst et all 120061: 
ISturm et al.ll2006h . we find that the mid-IR/X-ray 
luminosity ratio is characterized by a wide disper- 
sion and does not depend on the AGN class or 
the X-ray luminosity. The only exception is in an 
trend observed for the SF class to have a decreas- 
ing mid-IR/X-ray luminosity ratio at increasing 
X-ray luminosities. Some models predict smaller 
dust covering factors, or lower mid-IR over X-ray 
luminosity ratios, in more luminous sources. How- 
ever, according to these models, this effect should 
be observed in all AGN classes and at extreme X- 
ray luminosities, but the mid-IR/X-ray luminos- 
ity ratio in AGNls and AGN2s is constant and 



low ratios are also observed in SFs at moderate 
X-ray luminosities (~ 10 43 ' 5 ergs -1 ) as shown in 
Figure [8j 

In order to further investigate whether there 
is a dependency with luminosity, and whether 
previous studies did not cover a broad enough 
range of parameters (X-ray luminosities, spec- 
tral types) to observe such a trend, we produce 
a similar diagram using data from the litera- 
ture. Instead of the integrated mid-IR luminos- 
ity we use the monochromatic 6 /im luminosity 
and the absorption-corrected hard (2-10 keV) X- 
ray luminosity because available for all sources 
in the literature. The mid-IR over X-ray ra- 
tios as a function of X-ray luminosit i es from vari- 
ous AGN samples dBorys et al.l 20051 Sturm et al 



20061 : iLutz et all [20041 IWeedman et al.ll2006h are 
shown in Figure [9] In addition to the literature 
data, we report the average trend for the three 
AGN classes studied here. In this case, the average 
values were derived using the monochromatic 6 /im 
luminosity and the absorption-corrected hard (2- 
10 keV) X-ray luminosity as for the literature sam- 
ple. The resulting trends for our sample are con- 
sistent with the results shown in Figure [5] The 
literature data show a large dispersion, in agree- 
ment with our results, but no dependency on AGN 
type or X-ray luminosit y is observe d . Th ere is 
a group of sources from lBorvs et al ( 2005h with 
unusual low ratios and moderate X-ray luminosi- 
ties. These sources are sub-millimetre detected 
AGNs whose optical-IR emission is mainly dom- 
inated by the host-galaxy, similarly to the SF 
class. The X-ray emission of these sources is faint, 
L2- lOfcey^-lO 42 - 5 - 44 -^^- 1 , and they are 
characterized by lo w or moderate accretion rates, 
dM/dt ~0.01-0.80 (| Alexander et al.ll2005ah . 

The comparison with the ratios from the liter- 
ature samples confirms that the mid-IR/X-ray lu- 
minosity ratio in AGNs is characterized by a broad 
dispersion and there is no dependency with the 
X-ray luminosity. At low and moderate X-ray lu- 
minosities some AGNs are characterized by large 
mid-IR/X-ray ratios because of contribution from 
the host galaxy to the mid-IR luminosity. In SFs, 
the AGN mid-IR emission is generally fainter rel- 
ative to the absorption-corrected X-ray luminos- 
ity compared to AGNls and AGN2s. The host 
galaxy emission into the mid-IR can increase the 
mid-IR/X-ray ratios in SFs and make them consis- 
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Fig. 8. — Ratio between the mid-IR and the 
absorption-corrected broad-band (0.5-10 keV) 
rest-frame X-ray luminosity as a function of 
the absorption-corrected broad-band (0.5-10 keV) 
rest-frame X-ray luminosity. Symbols as in Fig- 
ure [6l The average ratios in bins of intrinsic X-ray 
luminosity derived for A sample (full symbols) are 
shown as diamonds and connected by solid lines. 
The colors correspond to the three SED classes, 
AGN1: blue, AGN2: red, and SF: green. 
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Fig. 9. — Ratio between the 6 /im and 
the absorption-corrected hard-band (2-10 keV) 
rest-frame X-ray luminosity as a function of 
the absorption-corrected hard-band (2-10 keV) 
rest-frame X-ray luminosity. Symbols refer 
to sources in the literature (reverse blue tri- 
angles and re d triangles : Se yfert Is and 



Seyfert 2s from Lutz et aL (12004), magent a di 



amonds: type 2 QSQ2s dSturm et al.l l2006h. or- 
ange squares: absorbed AGNs (jWeedman et al 
20061 ). and green crosses: sub- millimetre detected 



AGNs dBorvs et al J 120051) 1. The average ratios in 
bins of X-ray luminosity derived for sample A are 
shown connected by solid lines as full blue circles, 
full reversed red triangles, and green triangles, for 
the three SED classes, AGN1, AGN2, and SF, re- 
spectively. 



16 



tent or even larger than those observed in AGNls 
and AGN2s, but in some cases, the mid-IR/X- 
ra y ratios can be_100 times lower, as observed 



Borvs et al.l (|2005l ) sample. However, it is still 



not clear whether the observed low ratios are due 
to absorption of the mid-IR AGN emission or are 
intrinsic. Lower intrinsic ratios might imply dif- 
ferent dust properties (geometry, distribution, op- 
tical depth). For example, a low mid-IR/X-ray 
ratio could be indicative of lack of hot dust, as in 
the case the dust were located at a relatively large 
dist ance from the central h eating source, e.g. 10 
pc ( Ballantvne et al.1 120061 ). In this case, a de- 
pendency of mid-IR/X-ray ratio with the X-ray 
luminosity would not be expected. 

7. SEDs as a function of luminosity 

In SectionlH we compared the normalized SEDs 
of the three classes. Here, we investigate whether 
the average SEDs vary with the AGN luminosity. 
We use the absorption-corrected broad band X- 
ray luminosity, £o°5-io keV ^° separate each class 
(SF, AGN2, and AGN1) in three groups, sources 



< 10 44 ergs -1 , 2) L« 



with 1) L° %_ w keV 

44- ' 

1 The luminosity values that define the 



= 1()44~44.64 wga -l t and 3) L oorr_ wkeV 



corr 
5-10 keV 

> 10 44 64 



ergs 

three groups correspond to the value that sepa- 
rates Seyferts and quasars, £o°5-io keV = 10 44 
ergs -1 , and Lq°£1 w keV = 10 44 64 ergs -1 cor- 
responds to the median luminosity of all sources 
with L«£: 10 keV >l0 44 ergs -1 . In group (1), 4 
AGNls, 29 AGN2s, and 12 SFs are included; 10 
AGNls, 20 AGN2s, and 8 SFs in group (2), and 24 
AGNls, 10 AGN2s, and 2 SFs in group (3). The 
average SEDs of each class and group are shown 
in Figure [TOJ On the left side we show the SEDs 
regrouped by X-ray luminosity and on the right 
side by class. Note that the average absorption- 
corrected X-ray luminosities of each class within a 
group can differ by up to a factor of 2. The Elvis's 
median QSO template is shown for comparison. 
The average absorption-corrected X-ray luminosi- 
ties are shown as dashed lines at 0.5-10 keV in the 
rest-frame. 

The optical- near-IR (0.4-2/im) luminosities of 
the three classes in the lowest X-ray luminosity 
group (1) are all consistent within the la disper- 
sion (panel a of Figure [T0|) . However, the SEDs 
are characterized by different shapes, flat for the 



AGNls consistent with AGN emission at all wave- 
lengths, and optically red and flat in the IR for 
the AGN2s, consistent with being dominated by 
emission from the host-galaxy or by a reddened 
AGN in the optical and from the AGN in the IR, 
and optically red and blue in the IR for the SFs, 
consistent with being dominated by emission from 
the host-galaxy at all wavelengths. In the X-rays, 
AGNls and AGN2s show soft spectra, consistent 
with no absorption, while SFs show a hard spec- 
trum. The sources in this class are more than 10 
times less luminous than Elvis's QSO median tem- 
plate at all wavelengths. 

In group (2) (panel c), the optical-IR luminosi- 
ties of the three classes show similar properties to 
those in group (1). However, the difference in lu- 
minosity at short (A <0.4/im) and long (A >2^m) 
wavelengths in the three classes increases. The X- 
ray spectra are all similar and consistent with the 
canonical slope (ax— 0-7). The difference in lumi- 
nosity with Elvis's QSO median template is lower 
than in group (1), but they are still from three to 
ten times less luminous. 

In group (3) (panel e), the trend observed in 
the two previous groups is confirmed, the optical- 
IR luminosities of the three classes separate fur- 
ther, with the AGNls being from ten to hun- 
dred times more luminous than the SFs. In the 
X-rays, a large difference is observed between 
the absorption-corrected luminosities and the ob- 
served values in AGN2s and SFs, indicating large 
absorption, while no indication of X-ray absorp- 
tion is seen in AGNls. The AGN1 luminosities 
are now consistent with Elvis's QSO median tem- 
plate in the optical and X-ray, and higher in the 
IR. 

The average SEDs in each class show also a 
change in the overall shape as a function of X-ray 
luminosity (see right panels in Figure [T0|) . 

AGNls (panel b) have overall similar infrared, 
optical and X-ray SEDs at all luminosities, with 
only a slightly redder SED at A > l//m and less 
soft X-ray spectra at higher X-ray luminosities. 
The optical-IR luminosities scale with the X-ray 
luminosities, both uncorrected and corrected for 
absorption. 

The AGN2s (panel d) show a clear change in 
the optical-IR SED as the X-ray luminosity in- 
creases, with the overall SED becoming redder. 
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Fig. 10. — Left panels: Average rest-frame mid-IR- X-ray SEDs in Log(^L„) vs Log(z^) of all AGNs with 
absorption-corrected broad-band X-ray luminosity lower than 10 44 ergs -1 (panel a), io 44-44 - 64 erg s^ 1 (panel 
c), and >10 44 64 ergs -1 (panel e) are shown as filled circles connected by solid lines. Uncertainties correspond 
to the mean absolute deviation. Note that the X-ray SEDs are not corrected for absorption. The average 
SEDs for AGNs classified as SFs is shown in green, for AGN2s in red and for AGNls in blue. The dashed lines 
represent the average absorption-corrected 0.5-10 keV rest-frame X-ray luminosities for each AGN class and 
luminosity group. Right panels: Average rest-frame mid-IR- X-ray SEDs in Log(^L 1/ ) vs Log(^) of all AGNls 
(panel b), AGN2s (panel d), and SFs (panel f) with absorption-corrected broad-band X-ray luminosity below 
10 44 ergs -1 (cyan), io 44 ~ 44 ' 64 ergs -1 (magenta), and >10 44 64 ergs" 1 (purple). The number of sources 
used to derive the avera ge SEDs is anno tated. The black solid curve represents the median template for 



optically-selected QSOs ([Elvis et al.lll994l) 
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The most plausible explanation for this change is 
an increasingly contribution from hot dust contin- 
uum emission associated with the AGN, e.g the 
putative torus, or less contribution from the host 
galaxy. In the lowest luminosity group, the host 
galaxy still dominates, but at large luminosities 
its emission is significantly lower than that of the 
AGN. It is difficult to characterize the host galaxy 
emission in the mid-IR because it can change dra- 
matically for an elliptical or a late spiral or a star- 
burst galaxy. In all cases the host galaxy light 
peaks in the near-IR, ~l-2 nm, where the AGN 
light has a minimum ([Sanders et al.l Il989h . At 
longer wavelengths the host galaxy emission can 
decrease or increase, according to the galaxy type, 
and the AGN contribution increases. Thus, the 
total spectrum goes from blue if the host domi- 
nates to flat if the two contributions are similar to 
red when the AGN dominates. We can expect to 
see this trend for AGNs of increasing luminosity 
if the host galaxy luminosity is fixed or does not 
scale with the AGN luminosity. The X-ray spectra 
of the AGN2s change from being soft to hard as 
the luminosity increases. It is interesting to note 
that the AGN mid-IR luminosity does not scale 
with the absorption-corrected luminosity. Indeed 
the mid-IR luminosity of the most X-ray luminous 
AGN2s is less than twice that of those with inter- 
mediate X-ray luminosity, while the X-ray lumi- 
nosity is almost 10 times higher. This suggests 
mid-IR absorption in more luminous sources or a 
lower dust covering factor. 

The SF-like AGNs (panel f) show a similar 
change in SED shapes as the AGN2s. In this class, 
the near-IR (~1 /im) luminosity remains approxi- 
mately the same in the three groups, however the 
absorption-corrected X-ray luminosity is about 40 
times larger in the most X-ray luminous SFs than 
in the less X-ray luminous ones. The decrease 
in mid-IR luminosity relative to the absorption- 
corrected X-ray luminosity with increasing X-ray 
luminosities discussed in Section [5] is also visible 
here. This result is even more remarkable con- 
sidering that the host galaxy contributes to the 
mid-IR luminosity. 

The decreasing mid-IR luminosity compared 
to the absorption corrected X-ray luminosity ob- 
served in the most luminous AGN2s and in the SFs 
is equivalent to the low mid-IR over X-ray ratios 
found at high X-ray energies in SFs discussed in 



Section [6] and shown in Figures [8] and [9] From the 
analysis in the previous section, we concluded that 
the deficiency of mid-IR light is not dependent 
on the AGN X-ray luminosity. However, when 
the AGN X-ray luminosity is low the host galaxy 
might be luminous enough to compensate the mid- 
IR deficiency, preventing low mid-IR over X-ray 
ratios to be observable in less luminous AGNs. 
Assuming that a molecular torus surrounds the 
nuclear source, the deficiency of mid-IR light rel- 
ative to the X-ray luminosity observed in SFs and 
some AGN2s can be explained by lack of hot dust, 
larger optical depths, or large angles of the line of 
sight from the torus axis (jGranato fc Danesdll994 : 



Granato et al.lll997l ) 



8. The origin of AGN IR colors 



Lacv et al.l (|2004l ) identified a region in IRAC 



color-space (3.6-8.0 fim) where most of the 
AGNs detected i n all four IRAC bands are dis- 
tributed (see also Barmb v et al.l 2006; IStern et al 



12004 lHatziminaoglou et al.ll2005l ). AGNs occupy 
preferentially this region because of their red and 
almost featureless IR SED. Only five out of 87 
sources in sample A with detections in all four 
IRAC bands are not wi t hin t he 'AGN wedge' 
as defined by lLacv et aL ( 2004 ) . For fifteen ad- 
ditional objects, only an upper limit in one or 
more IRAC bands is available. We can there- 
fore conclude that Lacy's criterion is satisfied by 
80% of our X-ray sources in sample A observed 
with IRAC (this fraction is obtained with the 
pessimistic assumption that all the sources with 
IRAC limits would fall outside Lacy's wedge) or 
by a fraction close to 95% (97/102, if we assume 
that all the sources with upper limits lie within 
the wedge). Thus, this criterion appears to be 
fairly complete for a relatively bright X-ray sam- 
ple as the XMDS sample. Using the properties 
(SED class, luminosity, and absorption) derived 
for the AGNs in the XMDS, we investigate how 
these affect the IRAC colors and make them so 
distinctive for AGNs. 

In the hypothesis that different AGN types 
and column densities correspond to tori seen at 
different angles with respect to the line of sight 
and that the mid-IR radiation is dominated by 
the torus emission, a dependence of the IR col- 
ors on the AGN class (type 1 vs type 2) and 
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Fig. 11 — IRAC colors, F(8.0 /xm)/F(4.5 /im) versus F(5.8 /im)/F(3.6 ^m), as in lLacv et all (|2004 for our 
sample. Full symbols represent sources in sample A and open symbols represent sources in sample B. Each 
figures shows sources with different absorption-corrected X-ray luminosity (a), mid-IR (3-20 /im) luminosity 
(b), SED class (c), and column density as estimated in the X-rays (d). Each hgure is divided in 4 panels with 
sources in 4 different sub-groups as annotated. Left pointing arrows represent 5tr upper limits at 5.8 /im, 
downward pointing arrows represent 5 a upper limits at 8.0 /xm. The region enclosed by the dashed line 
corresponds to the region identified bv lLacv et al.1 ( 20041 ) where the majority of IRAC-detected AGNs lies. 
The dotted lines define three sub-regions within the Lacy's region, labeled Rl, R2, and R3, where the same 
number of sources in our sample lies. The number of sources in each group and the fraction in each region 
are listed in Table [2] 
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on the absorption in the X-rays is expected. 
In Figure [Til we reproduce Lacy's color dia- 
gram plotting separately sources in different in- 
tervals of absorption-corrected X-ray luminosity 
(panel a), mid-IR (3-20 /im) luminosity (panel 
b), SED classification (panel c), and column den- 
sity (panel d). IRAC colors are constrained for 
107 out of the 117 sources in sample A, and 
102 satisfy Lacy's color criteria: all AGNls, and 
most AGN2s (96%) and SFs (82%) We di- 
vided the area defined by lLacv et aL ( 2004 ) in 
three sub-regions, Rl, R2, and R3, that con- 
tain sources with increasing flux ratios or red- 
ness and, approximately, the same number of 
sources (~34). The equation of the line that di- 
vides regions Rl and R2 is Log(F(8.0)/F(4.5)) = 
0.32 - 1.25 x Log(F(5.8)/F(3.6)), and the equa- 
tion of the line that divides regions R2 and 
R3 is Log(F(8.0)/F(4.5)) = 0.661 - 1.25 x 
Log(F(5.8)/F(3.6)). In order to investigate 
whether any AGN property (class, X-ray lumi- 
nosity, mid-IR luminosity, and X-ray absorption) 
is correlated with the mid-IR redness, we derived 
the fraction of sources with a specific property in 
each region. The number of sources per region, 
including upper limits or detections only, are re- 
ported in Table H 

The position of the sources in the diagram is 
strongly correlated with the SED class. The frac- 
tion of AGNls increases from Rl to R3, while most 
of the SFs are in Rl and most of the AGN2s in R2. 
The observed dependency is the opposite of what 
expected by AGN tori models where redder SEDs 
are predicted in more obscured sources. A likely 
explanation for this correlation is the increasing 
host-galaxy contribution from AGNls to SFs. A 
correlation with the IRAC colors is also observed 
as a function of X-ray and mid-IR luminosities, 
with the fraction of the most luminous sources in- 
creasing and that of the less luminous decreasing 
going from Rl and R3. This correlation is consis- 
tent with less host contribution in the IR in more 
luminous AGNs. The least correlated parameter 
with the IRAC colors is the column density mea- 
sured in the X- rays. However, al l sources in the A 
sample outside Lacy et al. ( 2004 ) region are X-ray 
absorbed. 



9. Radio luminosity vs X-ray and mid-IR 
luminosities 

In star-forming galaxies, the radio and far- 
infrared (F IR; 42.5-122.5 Mm) emission are strongly 
corre lated ( Dickey fc Salpete"nll984 ; de Jong et all 
19851 ). The majority of AGNs (radio-quiet AGNs), 
about 90%, depending on the sample selection, fol- 
low the same correlation as star-forming galaxies, 
and a minority, i.e. radio galaxies and radio- 
loud quasars, shows an excess of radio emis- 
sion compared to radio-quiet objects with similar 



FIR emission (|Sopp fc Alexanderill991t iRov et al 
1998). The FIR over radio flux ratio is thus a 



good diagnostic tool to identify AGN activity. 
In recent studies, this correlation has been used 
by replacing the FIR emission with the 24 ^m 
emission because more easily av ailable for faint 



high-z sources thanks to Spitzer ([A pplcton et al 



12004 Iffigdon et alj 120051: iDonlev et al.ll2006l ). al" 

though the correlation between the 24 /im and 
the FIR emission is charact e rized by a large dis - 



persion (jBrandl et all 120061 : iBarger et all 120071) . 



The radio continuum is also correlated with the 



X-ray emission in AGNs (e.g.; iBrinkmann et al 



120001 ; ISimpson et al. 2006). A correlation between 
the radio and the X-ray emission is also observed 
in star- forming galaxies, however the X-ray lu- 
minosity is on average 400 times l ower than in 



AGN s of similar radio luminosity (jRanalli et al 
20031 ). Here, we investigate whether our data are 



consistent with the well known FIR-radio and X- 
ray-radio correlations and whether radio data can 
help identifying heavily obscured AGNs. 

There are 32 sources in our sample detected 
at 1.4 GHz, six belong to the B sample (three 
with dubious optical counterparts and three with 
poorly sampled SEDs). Upper limits to the radio 
flux are available for 68 sources (2 are in the B 
sample). Thus, there are 26 radio-detected and 
66 non-detected sources in sample A. Radio data 
can be used to identify the presence of an AGN 
by looking for an excess of radio emission com- 
pared to the optical and/or the IR flux. Ac- 
cording to IStocke et~all (|l992j ). radio- loud AGNs 
are characterized by a radio over optical ratio, 
R = F 5GHz /F 25QQ ^ greater than 10C@. We de- 
rive the radio flux at 5 GHz assuming a spectral 



2 The definition of radio-loud based on the parameter 
R=FsgHz/F • depends on whether the parameter R 
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Fig. 12. — Ratio between the rest-frame mid-IR 
(3-20 /mi) luminosity and the rest-frame 1.4 GHz 
radio luminosity vs the rest- frame 1.4 GHz radio 
luminosity. Symbols as in Figure El The thick 
dashed line represents the L(mid-IR)-L (Radio) 
correlation derived from the L(FIR)-L(Radio) 
correlation observed fo r radio-quiet AGNs and 
star-forming galax ies ( Sopp fc Alexander! Il99lt 



Ranal li et al. I l2003h . The thin dashed lines repre- 



sent the dispersion associated with the conversion 
from the FIR luminosity to the mid-IR luminos- 
ity. The thick dotted-dashed line represents the 
L(mid-IR)-L(Radio) correlation derived from the 
L(FIR)-L (Ra dio) correlation observe d for radio- 
loud AGNs (|Sopp k Alexander! Il99lh . The thin 
dashed lines represent the dispersion associated 
with the conversion from the FIR luminosity to the 
mid-IR luminosity. Plus signs represent 9 sources 
classified as radio-loud based on the radio over op- 
tical or radio over infrared flux ratio, crosses are 
over plotted on the 3 sources that satisfy both cri- 
teria (see text). The upward-leftward arrows rep- 
resent the median radio luminosity and luminos- 
ity ratio obtained from the non-detected sources 
in each group. 
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Fig. 13. — Ratio between the absorption-corrected 
rest-frame broad (0.5-10 keV) X-ray luminosity 
and the rest-frame 1.4 GHz radio luminosity vs 
the rest-frame 1.4 GHz radio luminosity. Sym- 
bols as in Figure [T2j The dashed line repre- 
sents the L(X-ray)-L(Radio) correlation derived 
for r adio quiet AGNs in th e ROSAT-FIRST sur- 
vey (jBrinkmann et al"1l2000h . The dotted line rep- 
resents the L(X-ray)-L(Radi o) correlation derive d 
for star-forming galaxies (jRanalli et al. 2003). 
The black solid diagonal line corresponds to a fixed 
X-ray luminosity of 10 42 ergs -1 . 
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index a^=0.8 (F„ oc v~ ai< \ and the optical flux 
by interpolating the observed data at the wave- 
length corresponding to 2500A rest-frame. Seven 
sources have R >100 (XMDS seq: 18, 58, 111, 
120, 403, 747, and 1219) of which two (XMDS 
seq: 403, and 1219) are optically blank and five 
are in sample B (XMDS ID: 111, 120, 403, 747, 
and 1219). We also calculated the parameter q 
defined as LogiFi^m /F '20cm) , and i dentified a ra 



dio ex cess in all sources with q <0.0 (jDonlev et al 
2005). Five sources satisfy this criterion, XMDS 



seq: 403, 738, 747, 1219, and 1246. Three sources 
are selected by both criteria (XMDS seq 403, 747, 
and 1219), for two (XMDS seq 111, and 120) the 
q parameter can not be derived because there is 
no IR coverage. In summary, the presence of an 
excess of radio emission remains uncertain in six 
sources for which only one of the two criteria is 
satisfied or can be applied, XMDS seq: 18, 58, 
111, 120, 738, and 1246, and it is confirmed by 
both criteria in 3 sources, XMDS seq: 403, 747, 
and 1219. 

In Figures fT2landfl~3l we investigate the correla- 
tion between the radio luminosity and the mid-IR 
and X-ray luminosity by comparing the ratio be- 
tween the rest-frame integrated mid-IR (3-20 ^m) 
and absorption-corrected rest-frame broad-band 
(0.5-10 keV) X-ray luminosities and the rest-frame 
radio luminosity at 1.4 GHz with the rest-frame 
radio luminosity at 1.4 GHz for the 32 X-ray 
sources that are detected at radio wavelengths. 
We show the luminosity ratio to remove the red- 
shift effect. The 9 radio-loud objects are shown as 
black plus signs, and crosses are over-plotted on 
the 3 more secure radio-loud AGNs. The 6 sources 
in sample B are shown as open symbols. Al- 
though these sources are classified as radio-loud by 
the criteria described above, based on their radio- 
to-MIR and radio-to- X-ray luminosity ratios, it 
would be more appropriate to consider them as 
radio-intermediate AGNs. For clarity we do not 
report sources that are not radio-detected, but for 
each class we show the median upper limit to the 
radio luminosity and to the luminosity ratio ob- 
tained from the non-detected sources. 

In Figure I12| we test whether our sample 



1 Hz" 1 and 
two relationships are 



is k-corrected and on what spectral indexes are assumed 
in the optical and in the radio. Typical values are 
LogR <0,5 for radio-qui et and LogiJ >2 for radio-loud ob- 
jects llStocke et al.lll992t) . 



shows the same correlation between the far-IR 
(60 /xm) luminosity and the radio luminosity at 
5 GHz as ob served in local radio-quiet and radio- 
loud AGNs (|Sopp fc Alexander! Il99lh . We con- 
verted the monochromatic luminosity at 5 GHz 
to that at 1.4 GHz assuming a radio spectral in- 
dex ^ =0.5 (/ oc v~ aR ) as in ISopp &: Alexander! 
(|199ll ) and the monochromatic luminosity at 
60 j-im to the mid-IR luminosity using the follow- 
ing equation, £3-20^™ = (°- 4 ± 0.8) x i/L 60fim , 
derived from the templates that fit the SEDs 
of our sample. The expected relationships are: 
Log(L 1A GHz ) = Log{L 3 -20nm) - H-02 ± 0.8 
for the radio- loud sources and Log(L\AGHz) = 
Lo3(L 3 _20^m) — 14.42 ± 0.8 for the radio-quiet 
ones, where L1.4 GHz is in erg s 
^3-20 M m m ergs -1 . The 
shown in Figure [T2] as thick dashed (radio-quiet) 
and dot-dashed (radio-loud) lines. The thin lines 
represent the dispersion associated with the con- 
version from vL^Q^m to £3-20^™- Most of the 
radio-quiet sources in our sample follow the ex- 
pected correlation. The Spearman correlation 
rank between the radio and the mid-IR luminosi- 
ties is ps=0.84 and the probability is lxlO -07 . 
The lowest luminosity ratios are observed for the 
sources considered radio-loud based on the opti- 
cal/radio or the IR/radio flux ratios, however they 
all show higher ratios than expected for radio-loud 
sources (dot-dashed line in Figure IT2"]). 

We find a similarly strong correlation be- 
tween the rest- frame radio luminosity at 1.4 GHz 
and the absorption-corrected rest-frame broad- 
band (0.5-10 keV) X-ray luminosity (Spearman 
rank correlation coefficient ps=0-91 and proba- 
bility of its deviation from zero 1.9xl0 -10 ). We 
compared the measured rel ationship with that 
found for radio-quiet A GNs ( Simpson et al. .2006; 
Brinkmann et al. 2000h (dashed line). The ex- 
pected relationship has been derived from the 
correlation observed between the monochromatic 
luminosities at 5 GHz and at 2 keV of a sample of 
radio-quiet (Lo qR < 1) sources from th e ROSAT- 
FIRST sample (|Brinkmann et al.ll2000l ). The ra- 
dio luminosity has been converted to L1.4 GHz 
assuming a radio spectral index of 0.7 and the 
X-ray luminosity at 2 keV has been converted to 
£0.5-10 keV assuming a spectral ind ex ax=0-& as 
assumed in lBrinkmann et a l. (2000). The final re- 
lationship is Log(L 1A GHz ) = Lo3(L .5-io kev) - 
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13.91 (see also ISimpson et al.l 120061 ) and is shown 
as a dashed line in Figure [13] The source in 
our sample show luminosity ratios consistent with 
those expected based on this relationship. 

The agreement with correlations derived from 
samples selected in different energy bands, and 
at different flux limits, reinforces the validity of 
these correlations. Since the radio luminosity 
can be considered an isotropic indicator of the 
nuclear unabsorbed X-ray luminosity for radio- 
quiet AGNs, these correlations can be used to 
probe the unabsorbed X-ray luminosity of radio- 
quiet AGNs and to identify absorption effects 
in X-ray observed AGNs, similarly to the anal- 
ysis based on the [O III] A5007 luminosity (see 



1997 



Risa 



Mulchaev et al 



Turner et al. 
iti et al.lll999T 



1994 
1997; 



Maiolino et al 
Bassani et al.l ll999). 



Alonso-Herrero et al 
1998; 



The XMDS sample offers an opportunity for 
such an application. Source XMDS ID 842 does 
not follow any of the correlations between the 
X-ray and the radio luminosities expected for 
AGNs or star-forming galaxies. Because of its 
low X-ray luminosity, lower than 10 42 ergs -1 , 
the evidence for AGN activity from the X-rays 
in this sources it is only based on the hard- 
ness of its X-ray spectrum. However, its radio 
over X-ray luminosity ratio is ~ 10 times larger 
tha n what expected for normal star-forming galax- 
ies (|Ranalli et al.ll2003l ), supporting the assump- 
tion that this source is an AGN. This behaviour 
suggests that the X-ray luminosity might be un- 
derestimated by a factor of 40, a hypothesis that 
would indicate that the source is Compton-thick 
and the observed X-ray emission is only scat- 
tered light. This source will be observed with 
the 7-ray observatory, INTEGRAL, as part of a 2 



Msec survey of the XMM-LSS field (jVirani et al 
20051 ). We estimate the luminosity and expected 



flux in the energy range 20-100 keV, where INTE- 
GRAL is sensitive, to determine whether INTE- 
GRAL will be able to detect it. Assuming that the 
absorption-corrected rest-frame 0.5-10 keV lumi- 
nosity is 40 times higher, the derived luminosity is 
£20-ioofceV=2.9xl0 43 ergs -1 , which corresponds 
to an absorption corrected flux ^o-ioo/ceV — 
7 x 10~ 12 ergcm" 2 s -1 , and an observed flux 
F20~iookeV — 6.8 x 10~ 12 ergcm -2 s -1 assum- 
ing a column density of 10 24 cm -2 . So far, IN- 
TEGRAL has not detected AGNs at these flux 



levels ( Beckmann et alj 2006; B assani et al . 2006), 
however the long XMM-LSS exposure might be 
able to reveal the heavily obscured AGNs in source 
XMDS seq 842 (aka Arp 54). 

10. Discussion 

10.1. AGN identification through SED fit- 
ting 

In this work, we have characterized the SEDs 



viouslv found (see e.e. Alexander et al.1 


200lfja 1 _ 


2002; 


Fadda et al. 2002:|Perola et al.l2004 


ilSzokolv et al. 


2004; 


Franceschini et al.l2005l Polletta et al.ll2006). 



the multi-wavelength SEDs of X-ray selected 
AGNs show a wide variety of properties. We 
divide the SEDs in three broad classes, AGN1, 
AGN2, and SF. The fraction of AGNs that fall 
into the various SED classes likely depends on 
the sample selection, i.e. spectral window and 
flux limit. For example, an optically selected 
sample will have a larger fraction of AGNls 
and an IR-selected sample will be biased against 
AGNs in the SF class. A different X-ray flux 
limit or a selection in a different X-ray band can 
also alter the fraction of sources in the various 
classes. 



For example 
finds 39% AGNls, 23% AGN2s and 
in a sample of AGNs with broad X-ray fluxes 
Fo.5-$keV > 10" 15 ergcm -2 s -1 . If the classifi- 
cation of the sources fitted with a Mrk 231 tem- 
plate in lFranceschini et al. I (|2005h is modified from 
AGN1 to AGN2 as assumed here, the fractions be- 
come 15% of AGNls, 42% of AGN2s, and 36% of 
SFs. The AGN1 fraction is lower than ours (15 
vs 33%), the fraction of AGN2s is slightly lower 
than ours (42 vs 50%) and the fraction of SFs 
is twice ours (36 vs 17%). These differences are 
likely due to our higher X-ray flux limit and re- 
quirement of detection in the hard X-ray band 
(2-10 keV). There is indeed an indication for the 
fraction of AGNls to decrease and the fraction of 
SFs to increase at fainter X-ray fluxes (see Paper 
II for a more detailed analysis). The trend of hav- 
ing less AGNls and more SFs at lower X-ray flux 
levels has to be confirmed and better quantified 
using similarly selected samples with analogous 
classification covering a broad range of fluxes. 

Although hard X-ray samples, as the one used 
here, provide the least biased AGN samples, some 



Franceschini et al 



(120051) 

38% SFs 
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AGN, e.g. Compton-thick AGNs are still missed. 
However, we do not expect the SEDs of Compton- 
thick AGNs to fall in a separate class. Known 
Compton-thick AGNs, both in the local universe 
and at high redshifts, show SEDs that are con- 
sistent with the AGN2 (e.g. Mrk 231, NGC 1068) 
and SF classes (e.g. NGC 4945, NGC 6240). From 
the SED analysis of a sample of X- r ay se lected 
Compton-thick AGNs, IPolletta et aD (|2006l ) esti- 
mates that 60% of Compton-thick AGNs fall in 
the SF class, and 40% in the AGN2 class. Thus, 
through SED fitting we can identify a large frac- 
tion of AGNs, even Compton-thick AGNs. The 
SED fitting technique provides an efficient AGN 
identification tool. The main advantage is that it 
is based on data that can be more easily obtained 
than X-ray data and spectroscopic data. More- 
over, this technique identifies the same or an even 
higher fraction of AGNs than the spectroscopic 
classification. Based on the sample of 49 spectro- 
scopically identified XMDS sources, we estimate 
that the SED fitting technique is able to identify 
AGN activity in 8 out of 11 (73%) AGNs that are 
spectroscopically classified as star-forming galax- 
ies, in 7 out of 12 (58%) AGNs that are spectro- 
scopically classified as type 2 AGNs, and in 25 
out of 26 (96%) AGNs that are spectroscopically 
classified as type 1 AGNs. Overall, we identify 
AGN activity in 40 out of 49 (82%) hard X-ray se- 
lected AGNs with available spectroscopic data and 
in 83% (50% of AGN2s and 33% of AGNls) in the 
entire X-ray sample, while optical spectroscopy re- 
veals an AGN in 38 sources (78%). The SED fit- 
ting technique can be applied to sources that are 
too faint at optical wavelengths for optical spectro- 
scopic observations, or with emission lines that are 
diluted by the host galaxy light or affected by ob- 
scuration dMoran et al. 2002 ; Barger et al. 2001 



Rigbv et al.l 120061 ) . and to sources in the redshift 
desert. There is evidence that these source might 
represent a significant fraction of AGN samples. 
Thus, SED fitting represents a powerful and effi- 
cient tool to identify AGN activity in samples for 
which optical spectroscopic observations are un- 
available or impossible to perform. 

However, this technique has also some limi- 
tations. Almost half of type 2 AGNs and a 
third of star-forming galaxy-like AGNs might still 
be missed and not identified as AGNs. More- 
over, the SED classification can be misled by the 



more dominant host-galaxy light at optical-near- 
infrared wavelengths as in many Seyfert 1 galaxies 
and thus classify a type 1 AGN as AGN2 or SF. 
These limitations might be overcome if well sam- 
pled multi- wavelength SEDs of large AGN samples 
in various classes with spectroscopic data become 
available. With such data the AGN templates can 
be improved and expanded and the SED fitting 
technique can be refined. 

In the next section we discuss the nature of the 
AGN that are not identified by the SED technique 
and by optical spectroscopy. 

10.2. The nature of elusive AGNs 

There are several examples of AGNs in the lit- 
erature whose optical spectra or broad-band SEDs 
resemble those of normal galaxies and are, there- 
fore, not identified as AGNs unless probed by X- 
ray or radio data. We refer to the se sources as elu- 
sive AGNs ( Maiolino et al. 20031 ). These sources 



belong to the SF class and might have many prop- 
erties in common with known categories of elusive 
AGNs found in the literature, e.g the optically 



dull discovered in optica l samples (IMaiolino et al 



2003 



Elvis et all Il98lt iRigbv et all 120051) . the 



X-ray Bright Optically Normal Galaxies discov- 
ered in X-ray samples (XBONGs; IComastri et al 



200% the type 3 QSOs found in IR - selected sam 
pies (|Martmez-Sansig re et alJ[20oi iLeipski et al 



2007). The rich taxonomy of these categories re- 



flects mostly the wavelengths at which they were 
discovered, while the scenarios proposed to explain 
their nature are similar. 

The properties of these elusive AGNs are usu- 
ally explained by two scenarios: 1) heavy ob- 
scu ration by dust with a large covering fac- 

tor dMarconi et al.l200oHFabianl2o"ollDudlev fc Wvnn- Williams! 
1997T T and/or 2) a radiatively inefficien t accre- 



tion flow (RIAF) (jYuan fc Naravanll2004l ). Large 
dust covering factors obscure the optical emission 
from the inner regions, thus no narrow emission 
lines are observed. The inefficient accretion ex- 
plains the lack of optical and ultraviolet radia- 
tion from the accretion disk and of emission lines 
from the broad and narrow line regions. The ex- 



isten c e of unabsorbed X BONGs (jSevergnini et al 



2003; iPage et alj 120031 ) in the X-rays disfavors 
the former hypothesis, although it is likely that 
both scenarios might be valid. Based on IR and 
sub-millimetre observations of some of these elu- 
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sive A G Ns dBorvs et al.l 120051: [Alexander et al 



l2005al |bl: iMartmez-Sansigre et all 120061 ). it has 



been proposed that the obscuring dust is associ- 
ated with a starburst and that these AGNs might 
be in an early phase of their growth. These objects 
might have already formed the bulk of their stellar 
mass, but the central BH has not fully grown and 



they are accreting at a l ow rate ([Alexander et al 



l2005ri iBorvs et alJliooil ). This last scenario com- 
bines the two described above. In summary, all 
these classes of AGNs are either explained by 
heavy absorption associated with dust with a large 
covering factor and/or by an inefficiently accreting 
black hole. 

We examine in more detail the properties of 3 
SFs for which spectroscopic data are also avail- 
able and fail in identifying the AGN, XMDS seq 
439, 487, and 842. Our goal is to find out what 
they have in common that might explain their elu- 
siveness and whether we can identify their AGN 
activity by other means in absence of X-ray data. 
All three sources are fitted with spiral templates, 
however XMDS seq 487 shows a small excess of 
emission in the near-IR with respect to the tem- 
plate. Such an excess is likely due to AGN-heated 
hot dust. Indeed, a fit with an AGN1 template is 
one of its secondary solutions. All three sources 
are absorbed in the X-rays (Nh>2x10 22 cm -2 ) 
and their X-ray luminosities cover a broad range, 
from 10 41 - 5 ergs" 1 to 10 444 ergs" 1 . The MIR-to- 
radio luminosity ratios of two sources are consis- 
tent with those observed in radio-quiet AGNs and 
normal galaxies, however XMDS seq 439 shows a 
lower than expected ratio (~13. 1 , compared to the 
expected 14.4; see Figure IT2"]) . This low ratio indi- 
cates radio emission in excess, and thus strongly 
suggests AGN activity. In summary, by using the 
MIR-to-radio luminosity ratio and the secondary 
solutions of the SED fits, we can identify AGN ac- 
tivity in 2 of these 3 sources. The AGN nature of 
the third source, XMDS seq 842, is actually not 
confirmed (see Sections [3.31 and 19")) . 

Overall, we do not find any evidence for ineffi- 
cient accretion in SFs. Indeed they can be intrinsi- 
cally powerful X-ray sources and their X-ray/radio 
luminosity ratios satisfy the same correlations ob- 
served for AGNls and AGN2s. The results from 
this study on these sources favor obscuration as 
the cause of the observed low mid-IR/X-ray lu- 
minosity ratios and elusive SEDs. However, an 



analysis of a larger sample with a broad range of 
luminosities and a rich multi-wavelength data set 
and, especially, more discriminating model predic- 
tions for various dust distributions and orienta- 
tions, linked with the X-ray properties, are neces- 
sary to investigate whether our results are simply 
consistent with an orientation effect, or indicate 
higher dust covering factors or larger dust inner 
radii. 

10.3. The missing AGN population 

The cosmic X-ray background is still unre- 
solved at high energies (>10 keV). The most re- 
cent background synthesis models predict that 
the missing population is made of sources at 
intermediate redshifts (z ~1), with hard X-ray 
spectra (r ~1.4), and moderate luminosities 



I^IO 43 e rgs~ 1N l (IWorslev et alJ 1200,4 iGilli et al 
20061) . ISimpson et all (|2006l ) suggest that the 



missing AGN population would be detected at 
radio wavelengths with fluxes of a few tens of 
/iJy. The expected radio luminosity of an AGN 
with an X-ray luminosity of 10 43 ergs _1 derived 
using the correlation shown in Figure [T2] is about 
10 29 ergs _1 Hz -1 . The observed radio flux, as- 
suming a redshift of 1 as predicted for the bulk 
of this missing AGN population, would be ~ 2.2 
/xJy and the observed 0.5-10 keV X-ray flux would 
be ~ 2.35x 10~ 15 erg cm~ 2 s _1 if unobscured and 
reduced by 20%, 50%, or 80% if obscured by a 
column density of 10 22 , 10 23 , or 10 24 cm -2 , re- 
spectively. The predicted radio fluxes are lower 
than the limits reached by the deepest radio sur- 
veys. Indeed, the r.m.s. of the deepest VLA sur- 
veys i s 2.7 AtJy in the 0.4 deg 2 SWIRE / Chandra 



field (jPolletta et alJ 120061 ; Owen et al 
prep.), and 5.3 uJv in the 0.04 d eg 2 GOODS 



North field ((Morrison et all 120061 ). Thus, the 



faint radio population is not a good candidate 
for the missing AGN population. Similar con- 
clusions have been recently derived from stacking 
the X-r ay data of faint ra dio-selected sources in 
HDFN dBarger et al.ll2007h . 



11. Summary 

Using the large data set from the XMM-LSS, 
SWIRE, CFHTLS, and VVDS surveys in the 
XMDS field (Paper I), the properties (SEDs, lu- 
minosity and absorption) of a hard X-ray selected 



2G 



sample of 136 X-ray sources are investigated. Pho- 
tometric redshifts are estimated and spectral types 
are determined using the Hyper-z code and a li- 
brary of galaxy and AGN templates using optical 
and infrared data (Section [3]). In this work, we in- 
vestigate the role of the main parameters responsi- 
ble for the observed properties of AGN SEDs, host 
galaxy contribution, gas obscuration, dust absorp- 
tion, and AGN power. Our analysis is based on 
multi-wavelength data, on the characterization of 
the average SEDs, and on the comparison of the 
luminosities emitted in the main energy windows, 
optical, mid-IR, X-rays and radio. Out of 136 hard 
X-ray sources, 132 are AGNs, two are clusters and 
two are star-forming galaxies. From the AGN 
sample, we select a sub-sample of 117 AGNs with 
high quality SEDs. The main results obtained 
from the analysis of this sample are summarized 
below. 

• AGN SED classification We fit the 

optical-IR (from 0.37 to 24/im) SEDs of the 
whole sample (134 sources, excluding the 
two clusters) to estimate photometric red- 
shifts and classify the sources. Spectroscopic 
data are available for 33 sources. We obtain 
photometric redshifts with <r(l+z)=0.11 and 
6% of outliers. The AGN sample is divided 
in three classes based on their SEDs, type 
1 AGNs (AGNls) (33±?%), type 2 AGNs 
(AGN2) (50i^%) and star-forming galaxy- 
like AGNs (SF) (17±jj%) (Section EI]). 

• Average SEDs We derive average SEDs, 
from hard X-ray to mid-IR wavelengths, 
for each AGN class (Section @J. AGNls 
show AGN-dominated SEDs from optical 
to infrared wavelengths, and mostly unab- 
sorbed X-ray spectra. Their optical spec- 
tra show a broad range of extinction, up to 
Ay=0.55, indicating that the AGN1 class 
includes some reddened AGNs. In AGN2s, 
hot dust emission associated with the AGN 
is detected at mid-infrared wavelengths, and 
emission from the host galaxy dominates in 
the optical. They show a broad range of X- 
ray properties, with the majority being ab- 
sorbed in the X-rays. Sources classified as 
SFs do not show any AGN signature at op- 
tical and infrared wavelengths, they are all 
absorbed in the X-rays. 



• Average SEDs vs X-ray luminosity and 
SED class The change in the average SEDs 
per class as a function of X-ray luminosity 
is analyzed in Section [7J The comparison of 
the average SEDs of each class shows that 
the optical-IR SEDs in the AGN2 and SF 
classes and the IR SEDs in the AGN1 class 
redden as the luminosity increases. This 
reddening is consistent with a lower rela- 
tive contribution from the host galaxy to the 
near-IR luminosity at larger X-ray luminosi- 
ties. 

• X-ray vs IR luminosity The majority 
of sources with high X-ray luminosities are 
AGNls, while SFs dominate at lower lumi- 
nosities. No correlation is observed between 
the X-ray or mid-IR luminosity and the X- 
ray absorption. 

We compare the mid-infrared (3-20 /nm) 
over absorption-corrected X-ray luminosity 
ratio with the X-ray luminosity and with 
observed values from various AGN samples 
in the literature. The luminosity ratios are 
characterized by a wide dispersion as the 
sources in the literature and we do not find 
any dependency on the X-ray luminosity. 
However, we find sources with significantly 
lower ratios in the SF class. Even lower ra- 
tios are found for sub-millimeter detected 
AGNs. These low ratios are explained by ei- 
ther obscuration in the mid-IR or by lack of 
hot dust as expected if the dust is located at 
a relatively large distance, e.g. 10 pc, from 



the c entral heating source (jBallantvne et al 
20061) . However, other explanations cannot 



be ruled out with the present data. 

• AGN mid-IR colors The analysis of AGN 
red mid-infrared colors as measured by 
IRAC reveals that redness is strongly corre- 
lated with the AGN class, but also with the 
X-ray and mid-IR luminosities and does not 
depend on the X-ray absorption (Section [H]). 
The correlation of the mid-IR redness with 
the AGN class is likely due to a decreas- 
ing host-galaxy contribution from SFs to 
AGNls. 

• Radio luminosity vs X-ray and IR 
luminosities Radio luminosities are com- 
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pared with the X-ray and mid-infrared lu- 
minosities cojffirmingjDre^^ 

tionships (jBrinkmann et alfeOOOnSopp &: Alexander 



19911) for AGNs (Section |JJJ>. A sub-set of 



sources shows radio emission in excess com- 
pared with the optical and/or mid-IR flux, 
consistent with begin radio-intermediate. 
We propose the use of the radio luminos- 
ity to identify heavily obscured AGNs by 
estimating the unabsorbed X-ray luminosity 
in radio-quiet AGNs and by comparing the 
estimated flux with the observed X-ray flux. 

AGN identification through SED fit- 
ting The SED fitting technique presented 
here succeeds in identifying signatures in 
83% of the objects in the hard X-ray sam- 
ple. The remaining 17±l% of AGNs are 
fitted with star-forming galaxy templates 
(SF class). A comparison with the spec- 
troscopic classification available for a small 
sub-sample indicates that the AGN1 class is 
100% reliable and 62% complete. The AGN2 
class is 29% reliable and 58% complete, and 
the SF class is 33% reliable and 27% com- 
plete. The SED and the spectroscopic clas- 
sification agree in 53% of the cases. 

The nature of elusive AGNs AGNs in 
the SF class are elusive AGNs at optical 
and IR wavelengths, similarly to optically 
dull AGNs, XBONGs, and type 3 QSOs. 
The properties of these sources can be sim- 
ply explained by large obscuration in the 
mid-IR, however it is not clear whether the 
large obscuration is an effect of orientation 
or whether it implies larger covering factors 
and/or dust at larger distance from the nu- 
cleus. We do not find any evidence for inef- 
ficient accretion in these sources as they can 
be powerful X-ray sources. 

The missing AGN population We inves- 
tigated whether the AGNs still missing from 
X-ray surveys and responsible for the bulk 
of the CXRB at high energies (10-30 keV) 
can be detected at radio wavelengths. The 
expected radio flux, derived assuming an X- 
ray luminosity of 10 43 ergs _1 and the rela- 
tionship shown in Figure \13\ and z=l is ~ 
2.2 /ijy This flux density is lower than the 
limits reached by the deepest current radio 



surveys. Thus, they will still be missing even 
in radio surveys. 

The SED analysis of the AGN sample presented 
here indicates the existence of a large fraction 
of AGNs that would remain unidentified in op- 
tical and/or IR-selected samples (SFs). These 
elusive AGNs show a broad range of luminosi- 
ties and they are mostly X-ray absorbed. The 
evolution of this population needs to be con- 
strained and taken into account in all evolution- 
ary models. In backward evolutionary models, 
where SEDs are associated with different types 



of sources dXu et al.ll200ll 120031: ISilva et al.| [2004: 



Treister fc Urry||2005hT the inclusion of AGNs with 



SF-like SED will greatly affect the estimated con- 
tribution of AGNs at all wavelengths, especially 
in the mid-IR where AGN-dominated SEDs are 
much redder than SF SEDs. The identification 
of these objects would be important in estimating 
the contribution of accretion ene rgy to the mid-IR 
background and source counts (ISilva et all 12004 
iTreister et"alll2006t iBallantvne et al.ll2006h . 
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Table 1 

Luminosities, photometric z and classification of the hard X-ray selected sample 
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Sp. type 


10 22 cm~ 2 
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L m i<j-7R 
Lq 


Av 


SED type 
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<30.95 


44.59 


0.20 


N6090 


SF 


B 






1.683 






14.693 


44.91 


45.10 


44.85 


<31.10 


44.79 


0.50 


TQSOl 


AGN1 


B 


229 


XMDS J022406. 4-041830 


1.128 






54.869 


44.27 


44.69 


44.43 


<30.69 


44.54 


0.05 


N6090 


SF 





Table 1 — Continued 



scq 


Source Name 


Zphot 




Sp. type 


10 22 cm~ 2 




T corr 
BX 

erg s — 1 


T corr 
H X 


Ll.4GHs 
erg s^ 1 Hz -1 


L© 


A v 


SED type 


SED Class 


Comment 


232 


XMDS J022449. 2-041800 


0.581 






<0.086 


43.75 


43.75 


43.46 


<30.00 


44.01 


0.00 


Seyl8 


AGN2 




233 


XMDS J022456. 0-041725 


1.293 






1.114 


44.45 


44.50 


44.25 


<30.83 


44.83 


0.55 


TQSOl 


AGN1 


Bu 


242 


XMDS J022437. 8-041520 


1.051 






5.200 


44.07 


44.24 


43.99 


<30.61 


44.14 


0.55 


spil_4 


SF 




246 


XMDS J022415. 6-041416 


2.106 






2.382 


45.18 


45.22 


44.97 


<31.33 


45.84 


0.50 


QSOl 


AGN1 




253 


XMDS J022451. 9-041209 


1.686 






<0.341 


44.79 


44.78 


44.55 


31.35 


46.08 


0.15 


TQSOl 


AGN1 




255 


XMDS J022408. 4-041149 


1.929 






30.718 


44.70 


44.94 


44.69 


31.26 


45.10 


0.05 


Arp220 


SF 


B 






0.357 






4.156 


42.90 


43.15 


42.91 


<29.50 


43.16 


0.10 


Sd 


SF 


B 


258 


XMDS J022447. 4-041049 


2.628 






21.429 


45.08 


45.23 


44.98 


<31.55 


45.55 


0.00 


QSOl 


AGN1 




270 


XMDS J022449. 2-040841 


0.958 






2.446 


43.91 


44.03 


43.78 


<30.52 


44.27 


0.00 


Seyl8 


AGN2 




271 


XMDS J022509. 5-040836 


2.042 






<0.471 


44.99 


44.98 


44.74 


<31.30 


45.84 


0.40 


BQSOl 


AGN1 




272 


XMDS J022501. 6-040752 


0.797 






<0.120 


44.10 


44.11 


43.89 


<30.33 


44.98 


0.25 


BQSOl 


AGN1 




279 


XMDS J022421. 3-040607 


0.260 






1.172 


42.77 


42.94 


42.69 


<29.19 


44.00 


0.35 


QS02 


AGN2 




280 


XMDS J022417. 9-040606 


1.633 






16.445 


44.59 


44.80 


44.55 


<31.07 


44.65 


0.20 


TQSOl 


AGN1 




281 


XMDS J022503. 2-040538 


0.930 






0.061 


44.61 


44.61 


44.36 


<30.49 


44.57 


0.00 


Seyl8 


AGN2 




282 


XMDS J022452. 1-040518 


0.189 






<0.041 


42.95 


42.96 


42.66 


<28.88 


43.69 


0.20 


QSOl 


AGN1 




288 


XMDS J022421. 2-040351 


0.566 






1.298 


43.61 


43.73 


43.48 


<29.97 


43.36 


0.35 


spi2a_12 


SF 




291 


XMDS J022452. 0-040258 


0.269 






54.378 


42.71 


43.40 


43.14 


29.69 


43.90 


0.00 


QS02 


AGN2 


M 


330 


XMDS J022333. 0-041525 


2.271 






71.370 


44.73 


45.03 


44.78 




45.60 


0.00 


QS02 


AGN2 




351 


XMDS J022356. 5-041105 


0.927 






0.679 


44.06 


44.10 


43.86 




44.25 


0.00 


Seyl8 


AGN2 




359 


XMDS J022325. 3-040922 


1.154 






17.337 


43.93 


44.20 


43.95 




44.26 


0.00 


spil.4 


SF 




403 


XMDS J022742. 1-043607 


3.493 






32.273 


45.46 


45.60 


45.35 


32.38 


45.39 


0.20 


N6090 


SF 


Bu 


406 


XMDS J022732. 7-043544 


0.713 






7.860 


44.01 


44.26 


44.02 


30.75 


44.49 


0.55 


TQSOl 


AGN1 




414 


XMDS J022726. 3-043327 


3.666 






10.703 


45.50 


45.56 


45.31 


<31.88 


46.34 


0.10 


BQSOl 


AGN1 




416 


XMDS J022812. 2-043230 


1.668 






<0.335 


44.94 


44.94 


44.73 




45.00 


0.05 


BQSOl 


AGN1 


No IR 


420 


XMDS J022729. 2-043225 


2.357 


2.290 b 




<0.577 


45.14 


45.13 


44.89 


<31.41 


46.20 


0.00 


QSOl 


AGN1 




427 


XMDS J022758. 6-043112 


0.859 






<0.131 


43.91 


43.91 


43.68 


<30.40 


44.80 


0.00 


119254 


AGN2 


No IR. 


430 


XMDS J022737. 1-043031 


0.760 






0.537 


43.93 


43.98 


43.74 


<30.28 


43.90 


0.40 


Seyl8 


AGN2 




438 


XMDS J022756. 3-042905 


0.375 






12.691 


42.76 


43.15 


42.89 


<29.55 


43.58 


0.00 


Arp220 


SF 


No IR 


439 


XMDS J022746. 0-042853 


1.988 


1.368 d 




13.195 


44.15 


44.36 


44.12 


30.93 


44.03 


0.20 


Sd 


SF 




440 


XMDS J022748. 8-042820 


2.574 






0.509 


45.43 


45.42 


45.22 


<31.53 


46.21 


0.00 


QSOl 


AGN1 


No IR 


449 


XMDS J022815. 2-042617 


2.504 






21.945 


45.03 


45.20 


44.95 




44.60 


0.40 


spil_4 


SF 


B, No IR 






2.438 






20.886 


45.00 


45.17 


44.92 




44.57 


0.00 


Seyl8 


AGN2 


B, No IR 


453 


XMDS J022802. 3— 042546 


0.568 






5.450 


43.47 


43.71 


43.48 




44.35 


0.30 


QS02 


AGN2 


No IR 


a 7n 
4 1 U 


YTVTTIQ TnOOQH/l K (1/11 SIS 


1.886 






<0.41 1 


45.01 


45.00 


A A 70 

44. i y 




45.41 


0.00 


119254 




rs , IN o 1 rv. 






0.418 






<0.065 


43.39 


43.40 


43.19 




42.51 


0.00 


spil_4 


SF 


B, No IR, M 






1.082 






<0.176 


44.41 


44.41 


44.20 




43.88 


0.00 


122491 


AGN2 


B, No IR 


487 


XMDS J022643. 6-043317 


0.489 


0.308 a 


3 


5.507 


42.95 


43.25 


43.00 


29.76 


44.29 


0.20 


spil.4 


SF 


M 


498 


XMDS J022629. 2-043057 


1.903 


2.031 a 


1 


<0.466 


45.41 


45.40 


45.20 


<31.29 


45.82 


0.00 


QSOl 


AGN1 




503 


XMDS J022649. 3-042920 


0.723 


0.634 d 


3 


4.310 


43.62 


43.83 


43.58 


<30.09 


44.38 


0.00 


QS02 


AGN2 




505 


XMDS J022649. 0-042745 


0.084 


0.327 a 


2 


<0.054 


43.07 


43.08 


42.86 


<29.42 


44.01 


0.00 


Seyl8 


AGN2 




521 


XMDS J022658. 8-042321 


1.754 


1.325 d 


3 


2.550 


44.84 


44.93 


44.74 


30.93 


44.89 


0.05 


Scy2 


AGN2 




523 


XMDS J022622. 1-042221 


1.586 


2.006 b 


1 


<0.456 


45.40 


45.39 


45.11 


<31.28 


45.96 


0.00 


BQSOl 


AGN1 




551 


XMDS J022342. 0-043533 


1.128 






1.723 


44.44 


44.52 


44.28 




45.12 


0.35 


QSOl 


AGN1 




561 


XMDS J022424. 1-043228 


1.678 






<0.338 


45.05 


45.05 


44.72 


<31.10 


45.79 


0.00 


BQSOl 


AGN1 





Table 1 — Continued 



scq 


Source Name 






Sp. type 


10 22 cm" 2 


L BX 


t corr 
BX 

erg s — 1 


t corr 
H X 


erg s — 1 Hz ~~ 1 


Lmid— IB. 

y 


A v 


SED type 


SED Class 


564 


XMDS J022350. 7-043157 


0.224 






<0.044 


42.71 


42.72 


42.51 




43.57 


0.00 


Seyl8 


AGN2 


565 


XMDS J022356. 8-043115 


1.051 






0.774 


44.26 


44.31 


44.06 




43.98 


0.55 


spil_4 


SF 


567 


XMDS J022432. 4-043036 


0.588 






<0.087 


43.85 


43.86 


43.59 


<30.01 


44.21 


0.10 


BQSOl 


AGN1 


571 


XMDS J022330. 2-043004 


2.404 


2.666 a 


1 


5.506 


45.40 


45.45 


45.21 




46.31 


0.40 


QSOl 


AGN1 


577 


XMDS J022438. 9-042705 


0.188 


0.252 c 


2 


<0.047 


43.69 


43.70 


43.39 


29.51 


44.31 


0.10 


Seyl8 


AGN2 


578 


XMDS J022350. 7-042703 


1.033 






<0.165 


44.22 


44.22 


43.88 




44.69 


0.25 


QSOl 


AGN1 


602 


XMDS J022351. 2-042054 


0.097 


0.181 a 


2 


2.836 


42.41 


42.67 


42.42 




43.64 


0.00 


M82 


SF 


626 


XMDS J022326. 0-043534 


1.149 






1.191 


44.63 


44.69 


44.45 




44.91 


0.00 


N6090 


SF 


708 


XMDS .T022605. 3-045803 


1.476 






<0.276 


44.75 


44.75 


44.50 


<30.97 


44.77 


0.00 


Seyl8 


AGN2 


709 


XMDS J022606. 7-045722 


0.661 






<0.098 


43.91 


43.91 


43.60 


<30.13 


44.26 


0.00 


Seyl8 


AGN2 


710 


XMDS J022627. 4-045710 


2.320 






0.931 


45.55 


45.56 


45.35 


32.11 


45.97 


0.50 


TQSOl 


AGN1 


718 


XMDS J022615. 1-045355 


0.896 






0.900 


44.08 


44.14 


43.94 


30.64 


44.23 


0.45 


spil_4 


SF 


720 


XMDS J022628. 9-045252 


2.027 






<0.465 


45.14 


45.13 


44.88 


<31.29 


45.40 


0.50 


BQSOl 


AGN1 


731 


XMDS J022554. 1-044921 


0.534 






2.107 


43.26 


43.42 


43.17 


<29.91 


43.93 


0.55 


QS02 


AGN2 


738 


XMDS J022556. 1-044724 


0.898 


1.010 a 


1 


<0.160 


44.49 


44.49 


44.25 


31.92 


44.66 


0.00 


Seyl8 


AGN2 


739 


XMDS J022617. 1-044724 


0.184 


0.140 a 


2 


2.820 


42.11 


42.38 


42.14 


29.32 


43.56 


0.05 


spil_4 


SF 


742 


XMDS J022514. 3-044659 


1.615 


1.924 a 


1 


<0.425 


45.29 


45.28 


45.04 


31.34 


46.23 


0.10 


BQSOl 


AGN1 


743 


XMDS J022625. 2-044647 


1.556 






2.256 


44.74 


44.80 


44.56 


<31.02 


45.25 


0.15 


Mrk231 


AGN2 


746 


XMDS J022512. 6-044633 


0.219 






5.213 


43.04 


43.36 


43.10 


<29.02 


43.02 


0.55 


Sb 


SF 


747 


XMDS J022640. 4-044606 


1.248 






23.531 


44.49 


44.79 


44.53 


31.61 


45.17 


0.00 


119254 


AGN2 






1.248 






23.531 


44.49 


44.79 


44.53 


31.61 


45.17 


0.00 


119254 


AGN2 


748 


XMDS .T022610. 9-044550 


2.886 






12.547 


45.32 


45.42 


45.17 


<31.65 


45.79 


0.40 


TQSOl 


AGN1 


755 


XMDS J022600. 1-044412 


0.730 






<0.108 


44.06 


44.06 


43.83 


<30.24 


44.59 


0.10 


Mrk231 


AGN2 


760 


XMDS J022531. 4-044210 


0.833 


1.228 d 


3 


150.000 f 


45.31 


46.21 


45.96 


<30.78 


44.41 


0.40 


Seyl8 


AGN2 


779 


XMDS J022321. 8-045740 


0.637 






16.471 


43.93 


44.28 


44.02 




45.26 


0.00 


Mrk231 


AGN2 


780 


XMDS J022332. 0-045740 


0.963 






<0.151 


44.31 


44.31 


44.07 




45.04 


0.00 


Seyl8 


AGN2 


782 


XMDS J022326. 3-045708 


0.839 


0.826 a 


1 


<0.125 


44.12 


44.12 


43.98 




44.72 


0.00 


Seyl8 


AGN2 


787 


XMDS J022317. 9-045527 


0.981 






9.800 


43.91 


44.14 


43.90 




43.69 


0.15 


Seyl8 


AGN2 


788 


XMDS J022353. 7-045510 


0.958 






0.085 


44.58 


44.58 


44.36 




44.68 


0.40 


QSOl 


AGN1 


789 


XMDS J022329. 1-045452 


0.646 


0.604 a 


1 


<0.089 


43.84 


43.84 


43.56 




44.44 


0.00 


Seyl8 


AGN2 


800 


XMDS J022403. 8-045120 


0.874 






<0.134 


44.13 


44.13 


43.88 


<30.42 


44.66 


0.55 


Seyl8 


AGN2 


801 


XMDS J022344. 4-045120 


0.959 






1.167 


44.09 


44.16 


43.92 




44.24 


0.00 


Seyl8 


AGN2 


807 


XMDS J022333. 0-044924 


2.039 


2.302 a 


1 


2.460 


45.08 


45.12 


44.87 




45.53 


0.00 


BQSOl 


AGN1 


817 


XMDS J022354. 5-044815 


2.428 


2.458 a 


1 


<0.657 


45.23 


45.21 


44.91 




46.67 


0.00 


QSOl 


AGN1 


820 


XMDS J022319. 4-044732 


0.640 


0.293 a 


2 


0.185 


43.18 


43.23 


43.03 




44.27 


0.00 


QS02 


AGN2 


825 


XMDS J022330. 6-044633 


1.735 






3.157 


45.05 


45.12 


44.84 




45.34 


0.10 


Mrk231 


AGN2 


828 


XMDS J022318. 8-044616 


0.673 






<0.100 


43.75 


43.75 


43.53 




44.28 


0.05 


QSOl 


AGN1 


840 


XMDS J022330. 9-044235 


2.297 






13.666 


45.05 


45.19 


44.94 




45.39 


0.55 


Mrk231 


AGN2 


842 


XMDS J022402. 4-044140 


0.010 


0.0433 e 


3 


1.761 


41.26 


41.49 


41.24 


29.28 


43.71 


0.30 


Sd 


SF 


844 


XMDS J022343. 2-044105 


1.493 






60.833 


44.46 


44.83 


44.57 




45.36 


0.05 


119254 


AGN2 






1.857 






86.702 


44.69 


45.07 


44.81 




45.32 


0.25 


122491 


AGN2 


846 


XMDS J022317. 2-044035 


0.765 


0.842 a 


1 


<0.128 


44.27 


44.27 


44.11 




45.61 


0.45 


BQSOl 


AGN1 


1197 


XMDS J022720. 2-045738 


1.116 






15.507 


44.16 


44.42 


44.18 


<30.68 


44.95 


0.25 


122491 


AGN2 


1199 


XMDS J022651. 6-045714 


0.290 


0.331 a 


2 


0.249 


43.54 


43.60 


43.37 


<29.43 


43.85 


0.00 


QS02 


AGN2 



Comment 



M 



Bu 



B 
B 



Table 1 — Continued 



XMDS seq 


Source Name 


Zphot 


z BpEC Sp. type 


ivref f 
N„ ' 

, „22 2 

10 22 cm 2 




erg s — 1 




Li.4GJf z 
erg s _1 Hz -1 


L 


Ay 


SED type 


SED Class 


Comment 


1201 


A.1V1UO — U4fcooUo 


0. 168 




0. 295 


42.28 


42.36 


42.12 


<28.77 


42.39 


0.00 


ocy 18 






1219 


XMDS J022701. 6-045158 


2.112 




25.148 


44.66 


44.86 


44.62 


32.45 


45.13 


0.55 


M82 


SF 


Bn 


1226 


XMDS J022711. 7-045038 


0.946 




0.215 


44.81 


44.82 


44.59 


<30.50 


44.93 


0.55 


BQSOl 


AGN1 




1227 


XMDS J022736. 8-045033 


0.445 




1.147 


43.32 


43.45 


43.26 


29.91 


44.19 


0.00 


Sey2 


AGN2 




1231 


XMDS J022731. 9-044957 


0.741 




<0.110 


43.62 


43.62 


43.39 


30.85 


44.18 


0.30 


Scyl8 


AGN2 




1236 


XMDS J022729. 0-044857 


1.513 




14.984 


44.31 


44.52 


44.28 


<30.99 


45.51 


0.40 


119254 


AGN2 




1246 


XMDS J022712. 8-044636 


1.446 




<0.267 


44.74 


44.74 


44.48 


32.85 


46.02 


0.00 


BQSOl 


AGN1 




1247 


XMDS J022633. 1-044637 


1.197 




3.040 


44.56 


44.68 


44.42 


<30.75 


44.91 


0.00 


QS02 


AGN2 




1248 


XMDS J022725. 4-044619 


0.034 




<0.029 


41.30 


41.32 


41.09 


<27.31 


42.60 


0.00 


Sa 


SF 


M 


1252 


XMDS J022716. 0-044539 


0.590 




<0.087 


44.07 


44.08 


43.82 


<30.02 


44.55 


0.50 


BQSOl 


AGN1 




1264 


XMDS J022751. 3-044251 


1.694 




<0.343 


44.95 


44.94 


44.79 


<31.11 


45.42 


0.10 


BQSOl 


AGN1 




1265 


XMDS J022712. 6-044221 


0.232 


0.205 a 2 


2.526 


42.94 


43.18 


42.93 


29.48 


43.95 


0.55 


spil_4 


SF 





Note. — XMDS scq: X-ray catalog sequence number (Paper I; Paper II); z p hot- photometric rcdshift; z spec : spectroscopic rcdshift; Sp. type: spectroscopic classification (1: type 1 AGN, 2: type 2 
AGN, 3: star-forming galaxy). N|/ , effective column density (Paper II). Upper limits are derived from the minimum value of 2.6xl0 20 cm" -2 corresponding to the Galactic value, when no additional 
absorption was required; I^bx - logarithm of the not corrected for absorption rest-frame broad-band (0.5-10 keV) luminosity; L^° x r : logarithm of the absorption-corrected rest-frame broad-band (0.5— 
10 keV) luminosity; L'fjx^ logarithm of the absorption-corrected rest-frame hard-band (2—10 keV) luminosity. All X-ray luminosities have been derived assuming an absorbed power-law spectrum with 
photon index T— 1.7 and column density N^^. Broad-band X-ray luminosities have been derived from the broad-band fluxes and hard-band X-ray luminosities have been derived from the hard-band 
fluxes (sec Paper II). I^i.aghz'- logarithm of the monochromatic radio luminosity at 1.4 GHz rest-frame; l^m i d — iR- logarithm of the integrated mid-IR luminosity in the 3—20 fim rest-frame wavelength 
^ range; Ay: optical extinction derived from the best-fit assuming the prescription reported in Calzctti ct al. (2000); SED type: best-fit template (see Figure^; Class: AGN spectral classification (SF, 

CTa AGN1, AGN2); Comment: B indicates a source in the B sample, Bu marks those in the B sample with unreliable poorly sampled SED. M refers to a source for which the absolute magnitude criterion 

for the photometric rcdshift determination was relaxed, No IR indicates sources without SWIRE coverage. 

Spectroscopic rcdshift from Garcet et al., in prep. 

b Spectroscopic redshift from,Gavignaud ct al. J20061) . 

c Spectroscopic redshift from Lacy ct al. (2006). 

d Spectroscopic rcdshift from fLe Fevre et all i2005|) . 

Spectroscopic redshift from NED. 

f Fixed maximum value of the intrinsic column density. 



Table 2 

Fraction of AGNs vs mid-IR color redness 



Group 


N 


% in Rl 


% in R2 


% in R3 


All 


107 (87) 


28 (25) 


32 (33) 


36 (36) 


L X < 10 1J - 8 


25 (24) 


37 (45) 


29 (34) 


11 (13) 


10 43 - 8 <L X < 10 44 - 3 


29 (21) 


37 (32) 


35 (34) 


16 (13) 


10 44 - 3 < L X < 10 44 - 9 


26 (20) 


23 (18) 


26 (28) 


26 (26) 


L X > 10 44 - 9 


22 (17) 


3 (5) 


9 (3) 


47 (48) 


imid-M < 10 1U b 


24 (17) 


33 (32) 


21 (21) 


18 (13) 


10 10 - 6 < L mld _ IR < 10 110 


29 (23) 


47 (50) 


35 (31) 


8 (10) 


10 110 < L mid _ IB < 10 112 


25 (23) 


17 (14) 


35 (41) 


21 (26) 


L mid -m > 10 11 2 


24 (19) 


3 (5) 


9 (7) 


53 (52) 


Normal 


(0) 


(0) 


(0) 


(0) 


SF 


14 (7) 


30 (27) 


9 (0) 


5 (3) 


AGN2 


55 (45) 


63 (64) 


62 (69) 


39 (35) 


AGN1 


33 (30) 


7 (9) 


29 (31) 


55 (61) 


N H < 10 2i ' 2u 


25 (24) 


33 (41) 


32 (38) 


11 (13) 


10 21 - 20 <N H < 10 22 06 


29 (25) 


33 (32) 


32 (38) 


21 (23) 


10 22 06 <N H < 10 22 - 83 


27 (16) 


30 (23) 


29 (21) 


21 (16) 


N H > 10 22 - 93 


21 (17) 


3 (5) 


6 (3) 


47 (48) 



Note. — The reported statistics refers to 107 sources for which IRAC obser- 
vations arc available. Of these 87 sources are detected in all 4 IRAC bands, 
for the remaining 20 sources upper limits were adopted, 4.3, 8.3, 58.5, and 65.7 
/j.Jy at 3.6, 4.5, 5.8, and 8.0 fim, respectively. The values in parenthesis re- 
fer to the sub-sample of 87 sources with detections in all 4 IRAC bands. Lx 
is the absorption-corrected luminosity in the 0.5—10 keV rest-frame interval in 
ergs -1 . L m id- 1 r is the mid-IR luminosity integrated between 3 and 20 fj.m 
rest-frame in Lq. NhIs the intrinsic column density density measured in the 
X-rays in cm -2 . 

Rl, R2, R3 represent three regions defined in color space by the flux ratios, 
F(5.8)/F(3.6) and F(8.0)/F(4.0). They arc shown in Figurc[TT] Only sources 
in the A sample arc taken into account. 
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